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I.  INTRODUCTION 


This  is  the  FY84  final  report  on  the  program  entitled  "Laser  Remote 
Sensing  of  Atmospheric  Pollutants"  supported  by  the  Air  Force  Engineering  and 
Services  Center  and  the  Air  Force  Aerospace  Medical  Division.  The  effort  is 
part  of  a  larger  ongoing  program  at  Lincoln  Laboratory  to  develop  laser  remote 
sensing  techniques  for  the  detection  of  chemical  species  in  the  atmosphere. 

The  specific  tasks  which  were  conducted  during  FY84  for  this  research 
program  consisted  of  the  following: 

(1)  The  continuation  of  the  CO2  DIAL  measurement  program, 

(2)  Further  development  of  the  Co:MgF2  DIAL  system, 

(3)  Initial  investigation  of  the  use  of  the  Co:MgF2  DIAL  system  for 
range  resolved  DIAL  measurements  of  atmospheric  species, 

(4)  Establishment  of  the  sensitivity  of  the  Co:MgF2  DIAL  system  for 
the  remote  sensing  of  atmospheric  HCl,  and 

(5)  Initial  investigation  of  frequency-shifting  techniques  for  the 
Co:MgF2  laser. 

Each  of  these  tasks  are  described  in  detail  in  the  following  sections. 
Supportive  documentation  is  included  in  the  appendices. 

II.  CO2  DIAL  INVESTIGATIONS 

The  study  of  the  use  of  CO2  differential -absorption  LIDAR  (DIAL)  for 
the  remote  sensing  of  atmospheric  pollutants  was  continued  during  FY84  and 
consisted  of  two  investigations.  The  specifics  of  these  two  investigations 
are  presented  in  Appendix  A  and  Appendix  B.  A  :ynopsis  of  the  studies  are 
presented  below. 


A.  Error  Reduction  in  Laser  Remote  Sensing:  Combined  Effects  of  Cross- 
Correlation  and  Signal  Averaging 

A  theoretical  analysis  was  made  of  the  extent  to  which  the  accuracy  of  a 
DIAL  measurement  may  be  improved  by  using  different  signal  processing  tech¬ 
niques,  in  particular  the  combined  effects  of  signal-averaging  and  temporal 
cross-correlation.  In  addition,  experimental  results  involving  signal  returns 
from  a  diffusely  reflecting  target  using  a  dual-C02  laser  DIAL  system  with 
both  heterodyne  and  direct  detection  were  found  to  be  in  excellent  agreement 
with  the  theoretical  analysis  (Appendix  A). 

The  DIAL  experiments  measured  the  effectiveness  of  using  the  ratio  of 
signal  returns  from  a  dual-laser  system  in  conjunction  with  signal  averaging 
in  improving  measurement  accuracy.  In  general,  significant  improvement  was 
obtained  by  first  averaging  the  individual  LIDAR  signals  and  then  taking  the 
ratios  of  the  averaged  signals  in  order  to  deduce  the  differential-absorption 
value.  This  is  in  contrast  to  taking  the  ratio  of  each  pulse-pair  DIAL  return 
and  then  averaging  this  ratio;  such  analysis  can  lead  to  a  significant  error 
due  to  the  exaggerated  effect  of  small  fluctuations  in  calculating  the  ratio 
value. 

The  degree  of  improvement  was  found  to  be  strongly  dependent  upon  the 
relative  values  of  the  standard  deviation  (fluctuation  level)  and  the  temporal 
cross-correlation  and  autocorrelation  of  the  LIDAR  signals.  However,  in 
general,  the  accuracy  of  the  DIAL  measurements  was  of  the  order  of  3  to  5%  for 
direct  detection  of  hard  target  returns  and  of  the  order  of  5  to  15%  for 
hetrodyne  detection.  It  should  be  added  that  these  results  are  primarily 
applicable  to  our  10-Hz  PRF  CO2  DIAL  system  and  should  not  be  construed  to 


apply  to  a  higher  PRF  system.  It  is  anticipated  that  a  higher  PRF  system 
could,  in  principle,  obtain  a  better  measurement  accuracy. 

B.  CO2  Laser  Remote  Sensing  of  Atmospheric  Ammonia 

The  CO2  DIAL  system  was  used  to  study  the  feasibility  of  the  laser 
remote  sensing  of  atmospheric  ammonia  (Appendix  B).  This  study  encompassed 
the  determination  of  the  optimum  CO2  laser  lines  for  the  remote  sensing  of 
ammonia,  the  first  demonstration  of  the  detection  of  atmospheric  ammonia  using 
a  CO2  DIAL  system,  and  the  observation  of  an  interplay  of  the  concentration 
of  water  vapor  and  ammonia  under  ambient  conditions. 

With  CO2  LIDAR  returns  obtained  from  topographic  targets  at  a  range  of 
2.7  km,  the  experiments  demonstrated  the  ability  of  the  DIAL  system  to  remote¬ 
ly  sense  the  presence  of  NH3  with  a  path-averaged  detection  sensitivity  of  5 
parts  per  billion  (ppb)  at  this  range.  The  concentration  of  atmospheric 
ammonia  was  observed  to  vary  between  a  negligible  amount  (less  than  5  ppb)  to 
values  as  high  as  20  ppb.  In  addition,  an  apparent  negative  correlation 
between  the  measured  concentration  of  ammonia  and  the  separately  measured 
relative  humidity  was  observed;  this  negative  correlation  may  be  due  to  the 
solubility  of  ammonia  in  water  vapor  aerosols  in  the  atmosphere. 

III.  Co:MgF2  LASER  DIAL  SYSTEM  DEVELOPMENT 

The  Co:MgF2  laser  reported  in  the  FY83  Final  Report  was  modified  in 
order  to  increase  laser  output  energy  and  spectral  tunability.  In  addition, 
the  DIAL  system  was  changed  to  provide  for  computer  control  of  the  laser  wave¬ 
length  and  DIAL  data  acquisition.  Some  of  the  details  of  the  laser  system  are 
presented  in  the  FY83  Final  Report.  The  specific  technical  modifications  are 
described  in  the  following  section. 
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Figure  1  shows  a  schematic  of  the  Co:MgF2  DIAL  system.  This  system  is 
similar  to  the  previous  system  with  the  following  changes:  (1)  the  solid, 
intracavity  etalon  is  now  0.25  mm  thick  and  coated  for  20%  reflectivity  at  1.6 
micrometers,  (2)  the  output  mirror  (M2}  was  changed  in  reflectivity  to  90% 
from  95%  to  reduce  the  internal  circulating  laser  energy  which  had  previously 
caused  optical  damage  in  mirror  Ml,  (3)  the  triple-element  birefrigent  tuning 
filter  (BTF)  was  angle  tuned  through  use  of  a  Burleigh  "Inchworm"  and  computer 
controlled,  and  (4)  a  60  cm  diameter  Cassegrain  telescope  was  used  instead  of 
the  off-axis  paraboloid. 

The  output  energy  of  the  Co:MgF2  laser  was  approximately  10  mJ  at  a  PRF 
of  3  Hz,  the  spectral  linewidth  was  0.15  cm"^,  and  the  pulse  length  was 
approximately  0.5  microseconds.  The  actual  laser  energy  transmitted  into  the 
atmosphere  was  estimated  to  be  about  7  mJ  per  pulse. 

The  temporal  characteristics  of  the  output  energy  are  shown  in  Fig.  2 
which  presents  a  photograph  of  the  oscilloscope  trace  of  the  free  running 
(not  Q-switched)  and  Q-switched  Co:MgF2  laser  output.  The  multi-spike 
output  under  free-running  conditions  is  typical  of  relaxation  oscillations  in 
a  solid-state  laser  under  nonequilibrium  pump  conditions. 

The  laser  energy  was  limited,  under  Q-switched  conditions,  to  less  than 
15  mJ/pulse  by  the  damage  threshold  of  the  back  mirror  (Ml  in  Fig.  1)  of  the 
laser  cavity.  It  is  anticipated  that  this  threshold  can  be  raised  by  an  order 
of  magnitude  by  changes  in  the  laser  cavity  and  laser  pumping  scheme;  these 
changes  are  currently  being  made. 

The  wavelength  of  the  laser  was  monitored  by  use  of  a  1-m  grating 
spectrometer  coupled  to  a  scanning  64-element  pyroelectric  detector  array 
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Figure  1.  Schematic  of  Co:MgF2  LIDAR  system. 
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Figure  2.  Oscilloscope  trace  of  Co:MgF2  laser  pulse. 
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placed  at  the  output  focal  plane  of  the  spectrometer.  The  array  provided  for 
wavelength  calibration  of  the  laser  on  a  single-shot  basis  with  an  accuracy  of 
approximately  0.5  cm"^.  Figure  3  shows  a  photograph  of  the  output  of  the 
spectral  array  detector  .  This  system  was  used  for  coarse  adjustment  of  the 
laser  frequency  and  to  ensure  stable  operation  of  the  Co:MgF2  laser.  Finer 
spectral  resolution  was  provided  by  a  fixed,  0.5  cm  spacing,  air-gap  Fabry- 
Perot  etalon  whose  output  was  recorded  by  the  computer  as  the  laser  frequency 
was  scanned;  this  provided  precise  calibration  of  the  wavelength  of  the  laser 
as  it  was  scanned. 

In  addition,  a  1-m-long  laboratory  absorption  cell  was  used  for  absorp¬ 
tion  calibration.  The  detected  normalized  absorption  signal  was  collected  by 
the  computer  system  for  comparison  with  the  LIDAR  (DIAL)  absorption  data. 

The  computer  was  used  to  control  the  wavelength  of  the  CorMgF^  laser 
and  to  collect  and  analyze  the  LIDAR  absorption  cell,  Fabry-Perot,  and 
normalization  (laser  energy /pulse)  signals.  Two  different  computer  programs 
were  written  which  were  used  for  path-averaged  and  range-resolved  DIAL 
measurements. 

For  path-averaged  measurements,  the  computer  scanned  the  CoiMgF^  laser 
in  small  increments  over  a  3  to  5  cm”^  frequency  region  and  recorded  the 
LIDAR  and  absorption  cell  data.  The  Fabry-Perot  data  was  used  to  calibrate 
the  wavelength  position  of  the  LIDAR  and  absorption  cell  data.  For 
range-resolved  data,  the  computer  "toggled"  the  laser  output  between  the  two 
wavelengths  chosen  so  that  they  were  on-resonance  and  off-resonance  to  the 
absorption  line  of  the  DIAL  detected  species.  The  range-resolved  LIDAR  data 
was  collected  by  a  transient  digitizer  and  averaged  over  a  selected  number  of 
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NUMBER  OF  ARRAY  ELEMENTS:  64 

CENTER  TO  CENTER  ARRAY  SPACING:  100  /dm 


~  25  cm  ^ 


Figure  3.  Output  of  spectrometer  pyroelectric  array  indicating  spectral 
distribution  of  Co:MgF2  laser  pulse. 


a 


I41S49-I 


Kli:  I  Srhfiii.itH-  i>r  (iualCOj  laser  (lirtfrcntial  al»s«»rplinii  lidar  system  with  capability  of  simultaneous  heterodyne  and  direct  detection 

of  backscattered  laser  sijjnals. 


hcttTddyne  and  direct  detection  of  the  same  hackscat- 
tered  lidar  signals, 

'I’he  1  )iAI ,  system  was  used  to  measure  and  record  the 
lidar  returns  from  a  diffusely  reflecting  target  using  both 
the  direct-  atul  heterodyne-detection  systems.  The 
target  was  located  at  a  2.7-km  range  from  the  labora¬ 
tory.  In  each  ex[)eriment,  a  total  of  6144  normalized 
f)ackscattered  lidar  signals  from  each  laser  was  recorded 
tor  later  stat  ist  ical  analysis.  'Phe  process  took  10  min, 
corresponding  to  a  pidse  re|)etition  frequency  (prf)  of 
-10  Hz. 

.■\  statistical  analysis  of  the  data  was  performed  and 
UK  hided  determination  ot  the  statistical  parameters 
presented  in  the  previous  section.  Specifically,  the  data 
were  iisefi  to  determine  ( 1 )  the  temporal  autocorrelation 
coetTicients.  p,,  and  (It)  the  temporal  cross  correla- 
t  ion  tor  |)ulse  pairs  ft,  .  for  longer  delays  p,,, ,  and  for  the 
signal-a\ eraged  pairs  /)„,  .  (H)  the  variance  of  the  indi- 
\idual  lasers,  rr  '  and  n\.  and  (4)  the  variance  of  the 
mean,  a;,,  and  a.y  From  these  parameters,  rr;;:  was 
computed  on  the  basis  of  Kqs.  (.7)  and  (H).  In  addition, 
after  averaging  the  signals  of  each  of  the  two  lasers  over 
n  pulses,  the  ratio  of  these  averages  permitted  a  direct 
measurement  of  rr,,  Comparison  of  the  computed  and 
measured  values  of  rr„,  yielded  a  direct  experimental 
deir’rminat ion  of  the  limits  ot  validity  of  Kqs.  (3)  and 
(.■>).  The  data  also  established  the  validity  of  Eq.  (8) 
anri.  most  important,  permitted  determination  of  the 
combined  eftect  iveness  of  signal  averaging  and  of  the 
dual  laser  cairabilitv  in  reducing  measurement  uncer¬ 
tainty  under  different  experimental  conditions.  Details 
ot  the  experimental  results  are  given  in  the  following 
sect  ion, 
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B.  Experimental  Results 

In  view  of  the  different  features  of  the  experimental 
results  obtained  under  various  experimental  conditions 
it  was  found  convenient  to  divide  the  direct  and  heter¬ 
odyne  DIAL  measurements  into  subgroups  dependent 
on  the  relative  numerical  values  of  tbe  standard  devia¬ 
tion  fT,(i  =  x,y)  and  the  cross-correlation  coefficient  p, . 
'I'he  subgroups  discussed  below  include  (1)  small  rr,. 
small  p, ;  (2)  large  c,,  large  p,.;  and  (3)  large  rr, .  small  p, . 
'I'his  division  was  chosen  because  the  different  DIAL 
measurements  indicate  it  is  these  relative  values  which 
determine  the  effect  of  signal  averaging  on  the  lidar 
returns.  A  fourth  subgroup,  namely,  small  ri,  and  large 
/», ,  is  not  considered  because  no  experimental  results 
obey  these  conditions;  this  is  associated  with  the  oh 
servation  that  uncorrelated  noise  often  dominate^  lidar 
fluctuations  at  low  rr,  levels,  so  that  a  large  /<,  \alue  i^ 
incompatible  with  a  small  n,  value. 

/  Small  Standard  Dri'latian  and  Snuill  (  ri,^^ 

(  Orndatinn 

Two  distinct  DI.AL  measurements  in\ol\ing  the  <ii 
reel  detection  ol  backstallered  lidar  return^  Imm  a 
diffusely  reflecting  target  had  small  n,  and  sm.ill  e 
values.  Under  these  experimental  conditiims,  it  ua^ 
generally  observed  that  rr,  and  rr,  ~  (i.'Jand  that  e  — 
0..1  0.4.  However,  while  the  ex|)erimental  condiiions 
were  the  same  in  the  two  experiments,  the  atmosphcrii 
condiiions  were  different.  In  one  of  these  experiments 
there  was  no  distinctive  variation  in  the  avtuage  value 
of  the  returns  over  the  10-min  experimental  pruiod. 
During  the  other  experiment,  there  was  a  slow  hut  al- 
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cross-correlation  coefficient.  'Phis  requires  deriving  the 
relationship  between  an  effective  signal-averaged 
temporal  cross  correlation  p,„  and  the  individual  pulse 
returns. 

Kquation  (.'!)  gives  the  approximate  relationship  be¬ 
tween  the  uncertainty  of  the  DIAL  measurement 
with  that  of  the  individual  laser  return  signals,  and 
IT, ,  and  the  pulse-pair  cross  correlation  p,..  The  terms 
in  Plq.  CD  refer  to  the  full  set  or  ensemble  of  data  points 
which  comprises  the  lidar  return  values.  If  one  assumes 
that  a  total  of  I'  data  points  have  been  measured  with 
each  of  the  la.sers,  a,  and  tt,  represented  the  normalized 
standard  deviations  of  the  1'  data  points  from  the  re¬ 
spective  lasers.  Signal  averaging  over  n  signal  return.-, 
where  n  «  1’.  then  yields  a  new  set  of  17ti  data  points, 
where  each  new  data  point  is  the  average  of  n  successive 
signal  returns.  One  may  then  formally  reexpress  Im). 
Cil  in  terms  of  these  signal-averaged  values  as 

where  a,,,  and  a,,,  represent  the  normalized  standard 
deviation  of  the  I'/d  signal-averaged  data  points  from 
lasers  x  and  respectively,  and  f>,„  is  the  cross  corre¬ 
lation  of  pairs  of  these  same  signal  averaged  data. 

It  is  particularly  important  to  note  that  since  tt,,,  and 
n,,,  have  been  determined  from  data  iiveraged  over  n 
pulses,  IT,,,  corresponds  to  the  normalized  standard 
deviation  of  the  ratio  of  the  returns  from  the  two  lasers 
after  averaging  these  lidar  signals  from  each  laser  over 
the  corresponding  T!  pidses.  'I’hal  is.  one  averages  t'irst 
and  then  takes  the  ratio.  This  averaging  then  acts  as 
a  low-pass  filter,  reducing  the  effect  ot  high-trecpiency 
fluctuations,  i.e.,  those  occurring  on  a  time  scale  that  is 
short  compared  with  the  time  recpiired  to  produce  the 
ri  pulses  being  averaged.  It  should  be  noted  that  using 
t  he  op[)osite  order  of  taking  t  he  rat  io  of  t  he  individual 
pulse  pairs  first  and  then  averaging,  while  permissible 
in  general.  Would  be  incorrect  in  the  context  of  Kcj.  (bt 
since  the  terms  on  the  right  han<l  side  ol  that  equal  ion 
refer  to  averaged  values.  In  addition,  taking  ratios  first 
may  seriouslv  bias  the  result-"'  (where  bias  refers  to  the 
flifferenie  between  the  tiu'an  value  of  t  he  est  itnator  and 
the  measured  value)  because  of  the  ai'centuated  effect 
of  small  fluctuations  of  the  individual  signals  on  (  heir 
ratio.  This  point  is  addressed  more  lull\  in  .Sec 
IV.H.J. 

The  terms  on  the  right-hand  side  of  Kig  (-'>)  tnav  be 
related  direitlv  to  the  individual  lidar  returns.  .Spe 
cificallv.  the  relationship  between  the  standard  devia 
lion  of  the  mean  ot  the  returns  from  lasers  I  and  a,,, 
and  n„  ,  and  their  res[)ecti\c  temporal  autocorrelation 
functions  is  given  b\  Kq.  (  1 1,  The  relal  ionship  ol  the 
effective  |)ulse  averaged  tenqioral  cross correlation 
(TH-fficient  I’,:.  .  with  the  pulse  [lair  value  /i,  and  with  the 
temporal  i  ross correlal ion  coefficients  ol  the  lasers  r  , 
tor  delay  t  imes  of  t  he  order  of  /  r  is  deriv  ed  in  .Appendix 
.-\  and  shown  to  be 


»l  .  ,  rt  . 


where  the  explicit  definition  of  p,,,  is  shown  in  Kq.  (An) 
to  be 


Ifsing  the  valuesof  a,,,  and  rr„  .  given  bv  Kq.  (-1),  one 
may  rewrite  Kq.  ((>)  as 


e,  -t-  'J  V  I  1  _  ;  ^ 


n  -  I 

1 

u  1 

1  +  -J  ^  (I  - 
-  =  1 

1  -f  ^  1 

' "  1 

1  -  I'n  1(1;  , 

Kquations  (fi)  and  (8)  relate  the  signal-averaged 
cross-correlation  coefficient  p„,  to  the  autocorrelation 
functions  of  the  individual  lidar  returns,  and 
and  to  the  cross-correlation  values  ft,  and 

Although  K(p  (8)  defines  the  functional  relationship 
between  p,„  and  the  fundamental  parameters 
I>,  ,  and  /<;n ),  it  should  be  noted  that  the  behavior  of 
these  fundamental  parameters  is  strongly  dependent 
on  experimental  conditions  and  difficult  to  predict  ac¬ 
curately  from  first  principles.  Measurements  which 
validate  K()s.  (b)  and  (8)  and  describe  the  behavior  of 
these  fundamental  parameters  in  different  experi¬ 
mental  conditions  are  |)resented  in  the  following  sec¬ 
tion. 

IV.  Experiments 

In  this  section,  we  present  the  results  of  experiments 
which  directly  measured  the  improvement  in  mea¬ 
surement  accuraev  t  hat  can  be  obtained  under  various 
experimental  conditimvs  by  vising  a  dvial -laser  DIAL 
system  in  conjunction  with  signal  averaging.  By  ex¬ 
perimentally  measuring  a,,,  and  /)„.  and  comparing  the 
results  obtained  with  those  predicted  by  the  theoretical 
analysis  of  t  he  jirevious  sect  ion.  we  established  the  va¬ 
lidity  of  K(|.  (8)  and  determined  the  limits  of  validity  of 
Kqs.  Ci)  and  t.b).  'I’he  analysis  also  prov  ides  a  theoret¬ 
ical  f  ramework  for  understanding  some  of  the  observed 
limitations  in  the  abilitv  of  signal  processing  to  reduce 
DIAL  measurement  uncertainties  in  these  experi¬ 
ments. 

A  Experimental  Apparatus  and  Technique 

.■\  schemat  ic  of  I  he  dual  ( '( )  •  laser  DIAL  system  used 
in  these  experiments  is  shown  in  Kig.  I.  .A  detailed 
de-scription  ot  (he  svstem  has  been  presented  previ- 
ouslv  '  '  The  svstem  emplovs  two  line-tunable  hybrid 
mini  I'K.-X  CO  lasers  '  w  hich  prov  ide  single-frequency 
pulsed  radiation.  I'he  two  lasers  are  separately 
triggered,  with  a  time  delav  At  between  tiring  of  the  two 
lasers.  T  he  value  ol  At  was  maintained  at  7  p sec  in 
the  experiments  desi  ribed  below;  over  this  short  time 
interval  the  atmosphere  is  effect ivt-lv  frozen.  The  beam 
paths  ol  the  two  lavers  are  joined  at  a  .atl/.AO  beam 
splitter,  rive  transmitted  beam  was  expanded  to  a 
b  cm  diameter,  and  the  receiver  is  a  2 -a  cm  Cassegrain 
telescope.  The  deteition  svstem  includes  two  inde¬ 
pendent  HgCd  I'e  detei  tor  units,  which  pi'rmit  both 
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employed,  much  shorter  decorrelation  times  have  been 
measured.  Experiments  carried  out  with  At  «  1  msec 
have  sl.own  that  the  value  of  p,  between  pulse-pair 
signal  returns  depends  strongly  on  both  the  nature  of 
the  reflective  target  and  the  detection  mode.-^  -'’  -'’  In 
general,  good  agreement  between  the  measured  value 
of  aj  and  that  predicted  by  Eq.  (3)  is  usually  obtained. 
However,  it  should  be  noted  that  the  derivation  of  Eq. 
(3)  involves  a  Taylor  series  expansion  and  is  correct  only 
to  first  order.  The  validity  of  the  equation  is,  therefore, 
limited  to  cases  where  higher-order  terms  are  negli¬ 
gible. 

B.  Signal  Averaging 

Signal  averaging  may  be  used  to  reduce  the  value  of 
the  standard  deviation  of  the  lidar  returns  from  each 
laser  in  Eq.  (3).  F'or  n  measurements,  the  normalized 
standard  deviation  of  the  mean  n,,  should  decrease  as 
assuming  the  measurements  are  independent. 
However,  experimental  DIAL  measurements  have 
shown  strong  departures  from  n“‘ behavior, indi¬ 
cating  that  successive  lidar  returns  were  not  indepen¬ 
dent  but  were  temporally  correlated.  The.se  experi- 
menUil  results  were  shown  to  agree  with  a  theoretical'-' 
study  which  established  the  relationship  between  n„ 
and  temporal  autocorrelation  as 

n  t  1  - 

(Tn  =  - H  1+2  ^  1  I  -  I 'n  ifi  •  (41 

\  n  [  .  =  t 

In  Eq.  (4),  n  is  the  normalized  standard  deviation  of  the 
full  set  of  indiv  idual  signals,  and  Py  is  the  tefnporal  au- 
t(K-orrelation  function  for  a  delay  time  jr.  where  t  is  the 
time  interval  between  pulses  of  each  individual  laser; 
the  autocorrelation  function  is  defined  in  Appendix  A. 
On  the  basis  of  Eq.  (4),  it  can  be  shown  that  small  tem¬ 
poral  autocorrelation  effects  among  successive  laser 
pulses  severely  limits  the  improvement  available  by 
signal  averaging."'  In  the  absence  of  temporal  auto¬ 
correlation.  F,q.  (41  reduces  to  the  usual  n"'-  be¬ 
havior. 

'I'o  avoid  confusion,  one  should  note*  the  difference 
between  the  cro.ss-correlation  coefficient  p,  in  Eq.  (3) 
and  the  autocorrelation  fund  ion  /),  in  E(i.  (4).  Cross 
correlation  refers  to  the  tem[)oral  correlation  of  the 
pulse-paired  returns  from  lasers  I  arid  '1.  while  auto¬ 
correlation  refers  to  the  correlation  ol  successive  pulse 
returns  from  each  individual  laser. 

C.  Physical  Processes  Relating  Lidar  and  Signal 
Processing  Parameters 

Before  continuing  with  the  signal  prix  essing  analysis, 
it  is  worthwhile  to  consider  the  pbvsical  mechanisms 
which  relate  the  lliict nations  (pf  the  return  lidar  signals 
/L.  /',  and  /' .  as  given  in  E<|s.  ( I )  and  (3),  with  the  signal 
prpxessing  parameters  of  E(ts.  (3)  and  (4).  Our  dis¬ 
cussion  will  be  ppialitative  rather  than  ipiantitative. 
tiuantitative  predii  tions  can  be  made  under  assumed 
ideal  londitions  but  are  beyond  the  scope  of  this 
paper. 

As  seen  from  E().  III.  the  statistical  and  temporal 
behavior  ot  the  return  signals  are  dependent  on  the 
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target  reflectivity  characteristics  p  and  the  atmospheric 
attenuation  effects  o,  each  of  which  has  associated 
statistical  fluctuations  and  temporal  correlation  prop¬ 
erties.  The  fluctuations  arise  as  the  result  of  speckle 
effects  at  the  target,  atmospheric  effects  in  the  form  of 
turbulence  or  relatively  long-term  atmospheric  ab¬ 
sorption  drifts,  and  exiterimentally  induced  fluctuations 
due  to  laser  beam  pointing  jitter,  laser  frequency  jumps, 
or  noise  processes  associated  with  the  optical  detection 
system. 

It  is  the  interplay  of  all  these  physical  effects  on  p  and 
tr  which  determines  the  statistical  and  temporal  char¬ 
acteristics  of  I’  and  F'  and  hence  the  values  of  ,  (Tv. 
p, ,  and  f)j.  Since  the  individual  effects  giving  rise  to 
fluctuations  will  have  different  spatial  and  temporal 
characteristics,  one  may  reasonaltly  predict  that  the 
relative  effectiveness  of  different  signal  processing 
techniques  will  depend  strongly  on  the  nature  of  the 
conditions  under  which  a  DIAL  measurement  is  carried 
out.  For  example,  heterodyne  detection  of  atmospheric 
lidar  backscatter  from  aerosols  results  in  short-term 
speckle-induced  fluctuations  of  the  reflectivity  term  p 
with  decorrelation  times  of  the  order  of  1  psec.--*  Under 
these  conditions  the  cross-correlation  coefficient  p,- 
would  be  expected  to  be  near  zero  and  should,  therefore, 
be  a  minor  factor  in  reducing  measurement  uncertainty. 
However,  during  the  time  period  required  for  signal 
averaging  over  a  large  number  of  pulses,  the  return 
signals  will  be  influenced  by  the  turbulence-induced 
fluctuations  in  o  with  decorrelation  times  of  the  order 
of  1  msec,----  ‘  or  the  even  longer-term  changes  in  o  due 
to  atmospheric  drift.  In  this  case,  p,  >  (I  and  can  lead 
to  a  significant  reduction  in  the  overall  uncertainty  of 
a  ratio  measurement  relative  to  the  uncertainties  of  the 
individual  measurements.  In  addition,  both  short-  and 
long-term  tem|)oral  effects  will  also  influence  the  au¬ 
tocorrelation  coefficient  /»,. 

Although  the  theories  of  speckle  and  propagation 
through  turbulence  are  well  understood,-'  it  remains 
difficult  to  take  all  the  factors  occurring  simultaneously 
into  account  with  sufficient  accuracy  to  make  an  accu¬ 
rate  prediction  of  these  values  o  priori  for  agiven  lidar 
system.  However,  these  theories  provide  a  basis  for  a 
qualitative  discussion  of  the  results  to  be  expected 
under  various  exiierimental  i-onditions.  The  lidar  ex¬ 
periments  de.st  ribed  in  this  paper  provide  information 
regarding  the  relative  effectiveness  of  signal  processing 
tei-hniques  lor  improving  the  accuracy  of  DIAL  mea¬ 
surements  under  different  experimental  conditions.  It 
is  also  hoped  that  the  results  will  provide  insights 
leading  to  a  more  complete  understanding  of  the  physics 
involved  in  laser  remote  sensing. 

III.  Analysis 

The  relationship  between  the  signal  return  ratios  and 
the  tempiiral  t  ross-corri'lation  (oi-lficient,  as  given  in 
K(|.  I  'll  lor  pulse  jiairs,  may  be  extended  to  incoriiorate 
signal  averaging,  (’ombining  both  temporal  cross 
correlation  and  signal  averaging  involves  expressing  the 
dual  laser  DI.-\L  relationship.  Kq.  (31,  in  terms  of  the 
signal  averaged  values  of  t  he  indiv  idual  lasers  and  their 
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the  theoretical  analysis.  These  results  are  used  to  study 
the  relative  merits  of  single-  and  dual-laser  lidar  systems 
and  to  help  establish  the  limits  of  accuracy  improve¬ 
ment  achievable  with  DIAL  measurements  by  the  use 
of  different  signal  processing  techniques. 

It  should  be  noted  that  this  theoretical  analysis  re¬ 
lates  quantities  which  are  experimentally  measured  and 
is  used  to  provide  insight  into  the  relative  magnitude 
of  the  effects  of  different  lidar  processes  on  the  overall 
accuracy  of  a  DIAL  measurement.  The  analysis  does 
not  predict  the  value  of  the  basic  quantities,  such  as 
lidar  fluctuation  levels  and  atmospheric  correlation 
times;  such  values  are  difficult  to  predict  from  first 
principles  except  under  selected  idealized  conditions. 

A  brief  review  of  our  previous  analysis  of  the  effects 
of  signal  averaging  and  temporal  cross  correlation 
considered  separately  is  given  in  Sec.  II.  A  theoretical 
analysis  which  combines  the  effects  of  signal  averaging 
and  temporal  cross  correlation  within  a  single  analvlical 
framework  is  presented  in  Sec.  III.  'Lhis  is  followed  in 
Sec.  IV  by  a  description  of  experiments  carried  out  with 
the  dual-CO-j  laser  DIAL  system  to  obtain  the  signal- 
averaged  variance  of  the  lidar  returns  from  the  indi¬ 
vidual  la.sers  and  of  their  ratio.  The  experimental  re¬ 
sults  are  shown  to  he  in  good  agreement  with  those 
predicted  on  the  basis  of  the  theoretical  analysis  given 
in  Sec.  III.  Results  are  also  used  to  help  determine  the 
limits  of  validity  of  the  theoretical  analysis  and  to  es¬ 
tablish  the  relative  merits  of  different  experimental 
techniques  for  the  various  experimental  conditions 
considered,  ('onclusions  ba.sed  on  the  results  of  the.se 
experiments  and  the  analysis  are  given  in  Sec.  V'. 

II.  Review 

The  lidar  equation  for  the  path-averaged  laser  remote 
sensing  of  a  molecular  constituent  in  the  atmosphere 
using  backscatter  from  a  hard  target  located  a  distance 
li  from  the  lidar  system  is  given  by 

/'■  =  I /' /\ //.t  ff/r ■  I  cxiil - +  (1  irt |.  (ii 

where  /',  is  the  instantaneous  received  optical  power  of 
l!:e  hackscattered  lidar  radiation  after  reflection  from 
the  target.  /’■  is  the  transmitted  power,  K  is  the  system 
optical  efficiency,  p  is  the  target  retlectivity,  A  is  area 
of  the  receiving  lelesi'ope,  n„  is  the  absorption  cross 
section  ot  the  miplei'ule  being  investigated.  is  the 
axerage  concentration  of  the  absorbing  molecule  over 
the  range  R .  and  o  is  the  atmospheric  extinction  coef- 
licient. 

Inal  )1,\L  measurement,  hackscattered  lidar  returns 
are  obtained  at  two  I requencies,  r  ami  c'.  whii  h  corre- 
s;>oMd  to  high  and  low-absorption  freijuencies  of  the 
molta  ular  spei'ies  investigated,  In  that  case,  the  con¬ 
centration  ol  the  spei  ies  is  dediu’ed  from  the  ratio  of  the 
lidar  returns  at  the  low-  and  high-abstirpt ion  trequen- 
(  ies  ol  the  species  investigated  rather  than  from  the 
indixidual  returns  according  to  the  relationship 

\  .  =  iiii  /’  /’I  _'i  .1  -  ,  i/i',  '  ji 

where  /’  =  I'r  I'l  is  the  received  signal  power  normalized 
to  the  transmitted  power,  and  the  primed  and  unprimed 


parameters  correspond  to  values  obtained  at  the  low- 
and  high-absorption  frequencies,  respectively.-' 

It  is  clear  that  fluctuations  in  the  lidar  signal  returns, 
P  and  P',  will  affect  the  accuracy  with  which  the  species 
concentration  may  be  determined  experimentally. 
It  was  noted  in  Sec.  1  and  the  experimental  data  will 
show  that  the  fluctuation  levels  of  P  and  P'  are  greatly 
influenced  by  the  detection  mode  employed.  This  is 
due  to  the  fact  that  in  direct  detection  the  return  lidar 
signal  is  proportional  to  the  square  of  the  electric  field 
amplitude  integrated  over  the  detector  field  of  view; 
that  is,  /’(direct)  «  /  E*EdA  \  in  this  ca.se  changes  in  the 
complex  phase  factor  are  eliminated.  However,  for 
heterodyne  detection,  P  is  proportional  to  the  integral 
of  the  product  of  the  electric  field  amplitude  of  the  re¬ 
turn  signal  beam  and  the  local  oscillator  field,  P(het- 
erodyne)  j'  Ei_,iE*dA.  'Fherefore,  phase  factors  have 
a  pronounced  effect  on  the  signal  fluctuation  levels  of 
heterodyne-detected  lidar  returns. 

Studies  have  been  made  of  the  extent  to  which  mea¬ 
surement  uncertainties  due  to  the  fluctuation  of  lidar 
returns  may  be  reduced  either  by  using  a  dual-laser 
system  to  take  advantage  of  cross-correlation  effects  or 
by  signal  averaging.  Results  of  the  studies  are  briefly 
reviewed  in  turn  in  the  following  section. 

A.  Temporal  Pulse-Pair  Cross  Correlation  in  a  Dual- 
Laser  System 

It  is  sometimes  possible  to  reduce  the  DIAL  mea¬ 
surement  uncertainty  by  taking  advantage  of  atmo¬ 
spherically  induced  temporal  cross  correlation  to  reduce 
the  fluctuations  of  the  ratio  of  these  returns.  However, 
to  achieve  this,  the  laser  beam  pulses  at  /'  and  v'  must 
be  separated  by  a  time  period  which  is  short  compared 
with  the  decorrelation  time  of  atmospheric  turbulence 
(~l  m.sec).--  'Phis  can  be  accomplished  by  using  two 
lasers  in  a  dual-laser  system  to  provide  two  nearly  si¬ 
multaneous  beams  oi)erating  along  the  same  path. 
With  lasers  I  and  2  operating  at  freipiencies  c  and  i'\ 
respectively,  it  has  been  found  that  the  normalized 
variance  of  the  ratio  of  the  return  signals  is  approxi- 
matelv  related  to  that  of  the  individual  return  values 
bv-’" 

,T  -  4t;  4  #T  ■  —  2  ft,  'T  ,  CT  .  (  M) 

where  n,  and  <t  are  the  standard  dex  iations  of  the  lidar 
return  signals  ol  lasers  1  and  2,  respectively,  normalized 
by  their  respect ixc  mean  x aloes,  a  is  the  normalized 
standard  deviation  ol  I  he  rat  io  ol  the  returns  ( t  =  .v/x  ). 
and  p,  is  the  temporal  pulse- pair  cross  correlat  ion  ol  the 
returns  Irom  the  txvo  lasers;  n  ,fT,,and  a  are  measures 
of  the  uncertainty  in  the  deterniinal  ion  nl  /'//*',  /’,  and 
/*',  respectively.  It  is  seen  in  Kq.  i:t)  that  a  positixe 
value  of  /I,  results  in  a  reduced  1)I.\L  measurement 
uncertainly  n,. 

Ivxperimenlallx  .  large  values  ol  /i,  haxc  been  obtained 
lor  pulse-paired  lidar  signals  rellecled  Iroin  slat  ionary 
I'argi-ts  xvhen  the  time  separation  Af  beixxcen  the  pulse 
pairs  is  ol  the  order  ol  iir  <1  msec.  .Atmospheric  (•Heels 
due  to  turbulence  max  be  considered  Irozen  oxer  this 
interval. '  When  a  moxing  target  '  or  source’’  is 
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Error  reduction  in  iaser  remote  sensing:  combined  effects  of 
cross  correlation  and  signal  averaging 


Norman  Menyuk,  Dennis  K.  Killinger,  and  Curtis  R.  Menyuk 


.A  sysienjrHir  unalysis  is  presenied  of  the  extont  lo  whiili  the  iieiurius  a  (iiHerent iiil-ahsnrpi ion  lidar 
(DIAL)  measurement  may  he  iniprt»ved  by  iisinj;  the eomhi ned  el t eels  of  siiiiial  it*. era^in^  and  temp«>ral  rross 
eorrelalioii.  Previous  studies  whieh  e«»nsidered  these  effeets  separalelv  are  extended  hv  inmrporatin^  both 
effeets  into  a  single  analytieal  framewt>rk.  In  adfiition.  experimental  results  itufilving  lidar  returns  from 
a  diffusely  relief  tin^'  tarj;et  usin^  a  laser  DIAL  system  v\iih  l>ofh  heterod> me  and  rjiret  t  deteetiim 

are  presented.  Lhese  results  are  shown  to  be  in  }i4n>d  agreement  with  t  he  theoretiejil  analysi^  and  help  estab¬ 
lish  the  limits  of  aeeuraev  aehieval>le  under  varit»us  experimental  i  tindii it»ns. 


I.  Introduction 

nit'ferential-atx'iorptioii  lidar  (Dl.AL)  systeni.s  have 
heeii  used  extensively  tor  the  single-ended  remote 
sensing  measurement  ot'  molecular  species  in  the  at¬ 
mosphere.'  Such  measurements  involve  use  of  at  mo- 
spheric  extinction  and  hackscatter  of  la.ser  radiation  at 
two  distinct  frequencies,  where  the  frequencies  are 
normally  chosen  to  coincide  with  high  and  low  absorp¬ 
tion  of  the  molecule  being  in\estigated.  I’he  concen¬ 
tration  of  the  molecular  species  is  deduced  from  the 
diflerence  in  the  measured  values  of  the  hackscattered 
laser  signal  returns  at  the  two  Irequencies. 

These  returns  are  subject  to  several  noise  sources, 
including  those  added  by  the  detection  process,  such  as 
shot  noise  in  the  case  of  heterodyne  detection  or  dark 
current  and  itackground  noise  for  direct  detection. 
These  noise  sources  generally  define  the  ultimate 
measurement  capability  of  a  DI.Al,  system.  ’  '  How¬ 
ever.  lidar  signals  traveling  through  the  atmosphere  are 
also  subject  to  significant  pulse-to-pulse  temporal 
tluctuations  due  to  several  factors  operating  on  the  laser 
beam  during  its  round  trip  from  transmitter  to  receiver. 
These  factors  include  extinction  and  .scintillation  effects 
ot  atmospheric  turbulence. '  '  speckle  or  glint  effects  at 
the  reflective  target.''*  and  laser  beam  wander  due  to 
either  atmospheric  turbulence  or  the  laser  source.'"  ' ' 


y  R.  Menviik  wilh  I  .1  Mar\l.tn<i.  I.atmraitirv  ti.r 

Plasma  &  Kuslim  Kner!;\  .Studio,  (  nllcac  Park.  Mar\land  _’oT  IJ;  ilu- 
other  authors  are  with  Mi  l'  I.im  oln  I.ahoratorv,  l.exuiKtoii.  Massa 
chusetls  (ViiTd-OoTS 
Received  S  Inly  IHS4 
0(H):t-69:i.=)  S.‘i  OKU  IS-U.SOJ  IKI  o 
c  198.0  Optical  Society  ol  ..Xmeru  a 

118  APPLIED  OPTICS  Vol  24.  No  1  Uanuary  1985 


In  addition,  since  the  heterodyne-  and  direct  detection 
modes  are  dependent  on  the  amplitude  and  phase  of  the 
lidar  signal  in  different  ways,  the  detection  mode  will 
also  influence  the  resultant  .SNR  and  fluctuation  level 
observed.-’  ^ In  many  cases,  notaltly  when  the 
average  signal  return  is  significantly  greater  than  the 
average  detector  noise  level,  it  is  the  pulse-to-pulse 
fluctuation  which  limits  the  measurement  accuracy  of 
the  average  lidar  return  signal  and  is,  therefore,  the 
factor  that  limits  the  accuracy  with  which  the  molecular 
species  concentration  can  he  determined  with  a  I)I.'\I. 
sy.stem.'"'" 

Previous  analyses  and  experimental  studies  investi¬ 
gated  the  extent  to  which  the  accuracy  ol  a  1)I.\1.  mea¬ 
surement  may  he  improved  through  use  of  I  wo  separate 
techniques,  namely,  by  signal  averaging  ol  the  lidar 
rettirns'"  '"  or  by  using  a  dual  laser  I)l.\l,  system  to 
exploit  the  fait  that  on  a  snllicientK  short  time  seale, 
both  laser  beams  may  eliectively  see  the  same  at  mo 
spheric  fluctuations.-"  In  this  paper,  we  extend  these 
studies  by  incorporating  both  techniques  within  a  single 
analytical  framework  which  effectively  combines  the 
influence  of  both  signal  averaging  and  temporal  cross 
correlation.  'I'his  extension  involves  a  theoretical 
analysis  which  relates  the  signal  averaged  statistical 
DIAL  measurement  error  and  the  temporal  autocorre¬ 
lation  and  cross  correlation  ol  the  individual  lidar 
pulses.  The  theoretieal  predietions  based  on  this 
analysis  were  compared  with  extensive  experimental 
mea.surements  obtained  using  both  heterodyne-  and 
direct -detection  COj  lidar  systems.  I'he  I fl.AL  exper¬ 
iments  measured  the  stat istical  error,  autocorrelation, 
and  cross  correlation  of  the  Dl.M,  returns.  'I'hese  ex¬ 
perimental  results  were  then  eom|)ared  with  those 
predicted  by  the  theoretical  analysis.  In  general,  the 
experimental  results  are  shown  to  agree  with  results  of 
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Appendix  A 


The  following  is  a  reprint  of  a  journal  article  published  in  Applied 
Optics,  1  January  1985,  entitled  "Error  Reduction  in  Laser  Remote  Sensing 
Combined  Effects  of  Cross  Correlation  and  Signal  Averaging". 
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nitrogen  which  produced  large  optical  distortion,  and  it  was  felt  that  a  more 
careful  design  of  the  cavity  would  alleviate  these  problems. 

VII.  SUMMARY  AND  CONCLUSION 

The  results  presented  in  this  report  have  described  progress  which  has 
been  made  in  the  development  and  measurement  of  the  1.7pm  Co:MgF2  DIAL 
system,  and  results  obtained  using  the  CO2  DIAL  system  for  the  remote 
sensing  of  atmospheric  ammonia. 

The  CO2  DIAL  measurements  indicated  that  NH3  may  be  present  in  the 
atmosphere  at  concentrations  as  high  as  30  ppb,  and  appears  to  exhibit  an 
inverse  dependence  on  the  water  vapor  aerosol  concentration.  The  investiga¬ 
tion  of  the  effects  of  signal  averaging  and  temporal  processes  in  the 
atmosphere  indicated  an  accuracy  on  the  order  of  a  few  percent  for  direct- 
detection  and  on  the  order  of  10  to  20  percent  for  heterodyne  detection. 

The  Co:MgF2  DIAL  system  has  been  shown  to  offer  significant  potential 
for  the  remote  sensing  of  HCl  in  the  atmosphere.  Further  investigations  are 
planned  using  calibration  sources  of  HCl  in  order  to  establish  the  sensitivity 
of  the  Co:MgF2  DIAL  system  for  both  path-averaged  and  range-resolved 
results. 
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returns,  one  deduces  that  the  path-averaged  sensitivity  of  the  DIAL  system  was 
approximately  0.2  ppm  of  HCl  over  the  3  km  range. 

VI.  INITIAL  INVESTIGATIONS  OF  FREQUENCY  SHIFTING  Co:MgF2  LASER 

Three  frequency-shifting  techniques  were  theoretically  studied  by  P. 
Moulton. 3  They  were  (1)  difference-frequency  generation  using  a  Nd:YAG  and 
Co:MgF2  laser,  (2)  Co:MgF2  laser  pumped  optical  parametric  oscillator 
using  a  AgGaSe2  crystal,  and  (3)  Raman  shifting  the  Co:MgF2  laser  to 
longer  wavelengths.  The  first  two  techniques  were  not  experimentally  investi¬ 
gated  because  (1)  the  temporal  overlap  of  the  NdrYAG  laser  pulse  and  Co:MgF2 
laser  pulse  is  very  small  since  pulsewidths  for  a  relatively  energetic 
Q-switched  Nd:YAG  laser  are  typically  10  ns,  which  is  much  shorter  than  the 
0-switched  Co:MgF2  laser  pulse  of  500  ns,  and  (2)  large  (>  0.5  cm)  AgGaSe2 
crystals  are  not  commercially  available  at  this  time.  Because  of  these 
limitations,  only  the  Raman-shifting  technique  was  experimentally 
investigated. 

Initial  Raman-shifting  was  attempted  using  the  Co:MgF2  laser  focused 
into  a  X/4  waveplate-i solated,  single-pass,  Raman  cell  containing  liquid 
nitrogen.  Results  of  this  experiment  were  negative,  because  of  the  low  power 
of  the  Co:MgF2  and  the  high-threshold  of  single-pass  stimulated  Raman 
scattering  (SRS). 

An  additional  experiment  was  attempted  to  observe  low  threshold  SRS  using 
liquid  nitrogen  in  an  optical  cavity.  Calculations  indicated  that  SRS  should 
have  been  above  threshold;  however,  no  signal  was  observed.  This  negative 
result  is  believed  to  have  been  caused  by  thermal  gradients  in  the  liquid 
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Figure  6.  Wavelength  scan  of  Co:MgF2  LIDAR  over  HCl  absorption  line 


Figure  5  shows  the  range-resolved  LIDAR  returns  as  the  laser  was  tuned 
on-  and  off-resonance  through  a  strong  water  vapor  line  near  1.7515  urn.  With 
an  expected  absorption  coefficient  of  8.3  x  10"^  (atm  cm)“^,  one  deduces 
from  the  observed  differential  signal  near  300  m  in  Fig.  5  that  approximately 
2.5  Torr  of  water  vapor  was  present.  This  number  is  in  approximate  agreement 
with  that  actually  measured  using  the  wet-bulb/dry-bulb  thermometer. 

Considerable  scatter  was  observed  in  the  range-resolved  returns  at  ranges 
greater  than  500  meters  due  to  the  low  signal -to-noise  ratio  for  the  absorbed 
on-resonance  signal.  Further  work  is  being  conducted  to  improve  this 
detection  range  by  increasing  the  power  of  the  laser  system. 

V.  REMOTE  SENSING  OF  HCl 

The  Co:MgF2  OIAL  system  was  used  for  the  path-averaged  remote  sensing 
of  HCl  in  the  atmosphere.  Laboratory  absorption  data  and  atmospheric  DIAL 
data  were  collected  as  the  Co:MgF2  laser  was  tuned  through  an  HCl  absorption 
line  near  1.7525  urn;  this  absorption  is  a  first  overtone  (v  =  0  +  2) 
transition. 

Figure  6  shows  the  absorption  spectrum  obtained  for  both  the  105-cm-long 
absorption  cell  containing  15  Torr  of  HCl  in  air  and,  simultaneously,  the 
LIDAR  returns  obtained  from  a  hillside  at  a  range  of  3  km.  One  deduces  from 
the  observed  60%  peak  absorption  due  to  HCl  in  Fig.  6a  that  the  absorption 
coefficient,  o,  for  this  line  is  0.45  (atm  cm)"^,  in  good  agreement  with 
previous  results. ^ 

As  evident  in  Figure  3,  no  obvious  HCl  absorption  feature  is  observed  in 
the  DIAL  spectrum.  With  a  5%  uncertainty  in  the  noise  background  of  the  DIAL 
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laser  pulses.  Each  LIDAR  return  was  normalized  to  the  transmitted  laser 
energy  per  pulse.  The  range-resolved  concentration  of  the  species  was  then 
computed  and  displayed  on  the  graphics  computer  terminal. 

IV.  INITIAL  RANGE-RESOLVED  Co:MgF2  MEASUREMENTS 

The  Co:MgF2  DIAL  system  was  used  to  measure  the  path-averaged  and 
range-resolved  concentration  of  water  vapor  in  the  atmosphere.  These  measure¬ 
ments  were  used  to  establish  the  capability  and  sensitivity  of  the  DIAL 
system. 

Path-averaged  DIAL  measurements  were  made  using  the  backscattered  DIAL 
returns  from  a  hillside  at  a  range  of  3  km  and  6.7  km.  Figure  4  shows  the 
DIAL  returns  as  the  laser  wavelength  was  tuned  through  several  water  vapor 
lines  near  1.767  um.  The  isolated  water  vapor  line  near  1,7679  urn  has  an 
expected  peak  absorption  coefficient,  a,  of  5.2  x  10"^  (atm-cm)"l  and  a 
linewidth  of  0,1  cm“1.2  the  observed  52%  absorption  of  this  line  in 

Fig.  4,  one  calculates  a  water  vapor  pressure  of  0.0023  atmospheres  (1.8 
Torn). This  value  is  a  factor  of  two  less  than  that  actually  present  as 
measured  by  a  wet-bulb/dry-bulb  thermometer  which  indicated  0.0048  atmospheres 
of  water  vapor.  This  discrepancy  in  the  determined  DIAL  H2O  concentration 
is  attributable  to  the  finite  linewidth  of  the  Co:MgF2  laser  (0.15  cm"M 
which  serves  to  reduce  the  peak  absorption  strength  and  to  broaden  the 
measured  absorption  line,  as  seen  in  Fig.  4.  Further  study  of  this  reduction 
is  planned. 

Preliminary  range-resolved  atmospheric  DIAL  measurements  were  made  with 
the  weak  backscattered  returns  from  naturally  occurring  aerosols  in  the 
atmosphere.  LIDAR  returns  were  detected  at  ranges  out  to  700  meters. 
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most  continuous  increase  in  the  magnitude  of  the  return 
signals  from  both  lasers  indicative  of  a  decreasing  at¬ 
mospheric  extinction.  The  results  obtained  under  both 
sets  of  atmospheric  conditions  are  described  in  the 
succeeding  sections.  These  experiments  were  per¬ 
formed  with  laser  1  operating  on  the  10.6-pm  P(22)  laser 
transition  and  laser  2  operating  on  the  lO.T-^m  P(28) 
transition. 

In  general,  it  is  desirable  to  have  the  frequencies  of 
the  two  lasers  as  close  together  as  possible  to  minimize 
the  decorrelation  due  to  turbulence-induced  scintilla¬ 
tion®'^  and  speckle  effects  at  the  target.  For  the  fre¬ 
quency  difference  between  the  P(22)  and  P(28)  tran¬ 
sition  and,  in  general,  for  other  frequency  differences 
within  the  range  of  CO2  lasers,  decorrelation  effects  due 
to  All  in  direct  detection  was  found  to  be  small  but  not 
negligible.  Specifically,  it  was  shown  that  Pc  ~  0.45  for 
i>  =  i'',22  while  for  i>  ^  v'  v/e  have  observed  Pc  *=  0.3-0.4. 
This  contrasts  strongly  with  decorrelation  effects  which 
occur  when  heterodyne  detection  is  used.  In  that  case, 
speckle  effects  at  the  target  cannot  be  averaged,  which 
results  in  a  strong  decorrelation  when  v  ^  v'.  This  is 
discussed  more  fully  in  Sec.  IV. B. 2. 

1.  Varying  atmospheric  extinction.  The  values 
measured  for  the  standard  deviation  of  the  mean  of  la¬ 
sers  1  and  2,  Onx  and  (T„y,  and  of  their  ratio  a„^  for  n  = 
1,2, 4,8, . . . ,  512  are  shown  in  Fig.  2(a)  for  the  case  where 
the  atmospheric  extinction  was  slowly  decreasing  over 
the  10-min  period  during  which  the  measurement  was 
carriedout.  The  values  of  ((r^j+  are  also  given 

since  they  represent  the  effective  standard  deviation  of 
the  measurement  in  the  absence  of  cross  correlation.  In 
effect,  these  values  represent  the  expected  standard 
deviation,  or  measurement  uncertainty,  of  a  DIAL 
system  which  uses  a  single  laser  switching  between  the 
on-  and  off-resonance  frequencies. 

There  are  several  features  of  interest  in  the  results 
shown.  They  are  as  follows:  (1)  the  reduction  of  Onx, 
(T„y,  and  a„(  with  n  is  significantly  slower  than 
(2)  for  n  =  1,  there  is  a  slight  improvement  in  mea¬ 
surement  accuracy  using  a  dual-laser  system  to  obtain 
the  ratio  value  (Tj  relative  to  (<Tj  -f  consistent  with 

a  small  positive  value  of  Pc  and  Eq.  (3);  and  (3)  the  im¬ 
provement  in  the  ratio  value  relative  to  (ffn,  -1 
is  seen  to  increase  continuously  with  increasing 
n,  indicating  that  Pn,  is  increasing  with  increasing  n. 

In  general,  the  measured  values  of  shown  in  Fig. 
2(a)  were  in  excellent  agreement  with  the  theoretical 
values  calculated  using  Eq.  (5).  In  addition,  the  func¬ 
tional  relationships  involving  the  autocorrelation  and 
cross-correlation  coefficients  were  found  to  agree  with 
the  theoretical  analysis  as  presented  in  Eqs.  (6)  and 
(8). 

Since  the  behavior  of  the  fundamental  parameters 
U’lt'  i>!\  -  l>i '  1  3''®  strongly  dependent  on  ex¬ 

perimental  conditions,  a  detailed  discussion  of  the 
physical  mechanisms  giving  rise  to  the  observed  func¬ 
tional  forms  of  these  parameters  is  considered  below. 

'I’he  slow  reduction  of  rr,,,  and  (t„,  with  increasing  n 
is  due  to  the  presence  of  temporal  autocorrelation,  Pjx 
and  /),, ,  in  the  lidar  returns  from  lasers  1  and  2,  re- 
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TIME  DELAY  III 


Fig.  2.  Lidar  signal  return  data  from  diffusely  reflecting  target  at 
a  2.7-km  range  using  direct  detection  during  period  of  continuously 
decreasing  atmospheric  extinction:  (a)  standard  deviation  of  returns 
from  lasers  I  and  2  and  their  ratio  as  functions  of  the  number  of  pulses 
averaged;  (b)  temporal  autocorrelation  and  cross  correlation  of  laser 
return  signals  as  functions  of  time  delay;  and  (c)  comparison  of 
measured  and  calculated  values  of  signal-averaged  cross-correlation 
coefficient  as  a  function  of  the  numlier  of  pulses  averaged. 


spectively,  as  shown  in  Fig.  2(b)  for  j  =  1, 2,4,8, . . . ,  512, 
where  the  temporal  cross  correlation  Pjxy  is  also  given. 
The  autocorrelation  function  is  unusually  large  in  this 
case,  apparently  due  to  the  presence  of  the  slow  atmo¬ 
spheric  drift  during  the  measurement  period.  This 
results  in  the  observed  long-term  temporal  correlation 
that  remains  essentially  unchanged  over  time  delays 
from  1 .6  to  50  sec.  The  measured  relationship  between 
the  standard  deviation  of  the  mean  for  the  two  sets  of 
la.ser  returns  and  the  autocorrelation  functions  is  in 
agreement  with  Eq.  (A6). 

As  noted  above,  the  increasing  ratio  of  relative  to 
(it);,  -t  iri;,)'/-  with  increasing  n  is  consistent  with  a 
value  of  p„r  that  increases  continuously  with  increasing 
n.  The  measured  values  of  p,„  vs  n  is  given  in  Fig.  2(c), 
and  just  such  an  increase  is  observed.  In  addition,  the 
values  of  p,„  were  calculated  on  the  basis  of  Eq.  (8)  using 
the  values  of  pjx,  Pjy,  and  pjxs  shown  in  Fig.  2(b)  and 
assuming  a  linear  intepolation  for  all  other  values  of  j. 
'I'he  calculated  values  of  p,„  are  also  given  in  Fig.  2(c) 
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and  are  seen  to  be  in  good  agreement  with  the  measured 
values.  The  reason  for  the  observed  increase  of  Pnc  with 
n  is  discussed  below. 

The  observed  variation  of  p„c  with  n  cannot  be  as¬ 
cribed  to  the  effective  low-pass  frequency-filtering  ef¬ 
fect  of  averaging  over  the  n  pulses,  since  that  effect  is 
similar  in  both  the  numerator  and  denominator  in  Eq. 
(8).  To  illustrate  this  point,  let  us  make  the  simplifying 
a  priori  assumption  pj^  =  Pjy.  Then  Eq.  (7)  becomes 

n-  1 

+  2  (1  -  j/n) II Its 

- (9> 

n  —  1 

1  +  2  ^{.l  -  j/n )pjx 
;  =  i 

The  experimental  results  of  Kjelaas  et  a/.®  indicate 
that  the  temporal  cross-correlation  coefficient  pjxy  is 
proportional  to  the  temporal  autocorrelation  coefficient 
pjx's  i.e.,  Pjxs  =  Kpjx.  However,  if  this  proportionality 
constant  is  extrapolated  back  to  a  zero  time  delay,  one 
would  obtain  K  -  p,.,  and  Eq.  (9)  would  then  become  Pnc 
=  p,.,  which  clearly  disagrees  with  experimental  re¬ 
sults. 

This  apparent  discrepancy,  which  arises  when  the 
proportionality  constant  is  extrapolated  back  to  zero, 
is  due  to  the  role  noise  plays  in  the  averaging  process  for 
autocorrelation  and  cross  correlation.  In  general,  lidar 
return  signals  contain  both  a  true  signal  and  a  noise 
component.  The  noise  components  of  the  two  lasers 
are  uncorrelated  with  each  other.  Therefore,  in  the 
averaging  proce.ss  given  in  Eqs.  ( A4)  and  ( A5)  the  effect 
of  the  noise  will  effectively  be  reduced  or  eliminated. 
The  effect  i)f  noise  will  he  similarly  reduced  or  elimi¬ 
nated  in  temporal  autocorrelation  measurements  by  the 
averaging  process  given  by  Flq.  ( A:{)  when  the  time  delay 
jr  is  greater  than  the  coherence  time  of  the  noise  com¬ 
ponent.  The  results  of  Ijelaas  vt  a/*’  indicate  that  this 
coherence  time  is  <2  msec;  for  delays  greater  than  this, 
the  expression  pjxs  -  Kl>ix  he  expected  to  hold  to 
within  experimental  uncertainties.  However,  for  zero 
time  delay,  i.e.Jr  =  0,  (/,(//,. -t-  jr))  in  Eq.  (All) 
becomes  (Ixiti;)’).  Since (0. ),  which  corresponds  to 
the  /eth  lidar  signal  return,  includes  a  noise  component 
Oidk)’)  will  yield  a  noise-squared  term  which  is  not 
eliminated  during  the  averaging  process.  Hence  the 
zero  time  delay  autocorrelation  coefficient,  which  is  by 
definition  equal  to  unity,  uniquely  contains  a  significant 
noi.se  component  after  signal  averaging.  'I'herefore,  p, 
<  Pjxs/ PjxsOT  K  =  j'Pi ,  where  C >  1.  With  this  further 
modification,  Eq.  (9)  becomes 
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Equation  (10)  gives  insight  into  the  physical  mecha 
nism  which  may  yield  a  rising  value  o(  p,„  with  in¬ 
creasing  n  as  long  as  py,,  >  0.  The  erpiation  indii'ales 
it  is  the  fact  of  H  multiplied  by  the  summation  term 
in  the  numerator  which  causes  this  increase.  I’herefore, 
a  summation  term  which  remains  positive  out  to  longer 


time-delay  periods  will  result  in  a  greater  increase  in  p„, . 
This  is  just  the  situation  which  occurs  in  the  case  of 
long-term  atmospheric  drift. 

A  necessary  condition  for  this  explanation  to  be  valid 
is  that  there  be  a  significant  increase  in  the  ratio  of  the 
cross  correlation  to  autocorrelation  on  going  from  a  zero 
time  delay  to )  =  1.  It  is  seen  in  Fig.  2(b)  that  there  is 
an  ~65-70%  decrease  in  the  autocorrelation  on  going 
from  zero  delay  to  the  delay  time  as.sociated  with  j  =  1 
(0.1  sec),  while  the  corresponding  decrease  for  the  cross 
correlation  is  only  80%.  This  is  in  accord  with  the 
model  proposed  (i.e.,  f  >  1)  to  explain  the  observed 
variation  of  p„,  with  n. 

It  should  be  noted  that  Eqs.  (9)  and  (10)  involved 
simplifying  assumptions  which  are  not  strictly  true.  In 
general,  as  we  have  seen,  Pjx  ^  Pjy,  and  there  is  con¬ 
siderable  scatter  in  the  experimental  values  of  Pjx ,  Pjy, 
and  pjxy  Therefore,  any  direct  comparison  with  ex¬ 
periment  must  be  made  using  Eq.  (8). 

2.  Constant  average  atmospheric  extinction.  A 
second  similar  experiment  carried  out  over  a  period 
during  which  the  average  atmospheric  extinction  re¬ 
mained  essentially  constant  led  to  the  results  in  Figs. 
3(a),  (b),  and  (c).  They  show  both  important  similari¬ 
ties  with  and  marked  differences  from  the  corre¬ 
sponding  results  described  above  and  shown  in  Figs. 
2(a),  (b),  and  (c),  respectively. 

Comparison  of  the  results  shown  in  Fig.  3(a)  with 
those  of  Fig.  2(a)  indicate  the  following  significant  dif¬ 
ferences  in  behavior:  (1)  although  the  reduction  of  Onx 
and  (T„s  with  n  remains  considerably  slower  than 
the  variation  is  significantly  more  rapid  than  was  ob¬ 
served  for  the  varying  extinction  case,  and,  (2)  the  im¬ 
provement  in  measurement  accuracy  obtainable  using 
relative  to  (a'ix  +  is  seen  to  increase  with 

increasing  n  to  n  =  8  and  then  to  decrease,  in  marked 
contrast  with  the  constantly  increasing  improvement 
with  increasing  n  observed  in  Fig.  2(a). 

The  more  rapid  decrease  of  rr,,.,  and  with  n  for  the 
constant  average  extinction  is  consistent  with  the  more 
rapid  decrease  of  the  autocorrelation  functions  with 
increasing  time  delay  j  for  this  case,  as  seen  in  Fig.  3(b). 
'I’he  relation.ship  between  the  standard  deviation  of  the 
mean  for  each  of  the  lasers,  a,,  ,  and  it,,,  ,  and  Py  ,  and  pjy, 
respectively,  are  in  agreement  with  Eq.  (A(i). 

The  measured  values  of  the  signal-averaged  cross- 
correlation  coeflicient  p,„  as  defined  by  Eq.  (A7)  is 
shown  as  a  function  of  n  in  Fig.  3(c).  In  addition,  the 
values  ol  /)„,  were  calculated  on  the  basis  of  Eq,  (7)  using 
the  values  of  p,,,  p,,,  and  py,,  given  in  Fig.  3(b)  and 
assuming  a  linear  interpolation  for  all  other  values  of 
J.  'I’hese  calculated  values  are  also  shown  in  Fig.  3(c) 
and  are  seen  to  he  in  good  general  agreement  with  the 
measured  values. 

The  occurrence  of  a  maximum  value  of  p,„  at  n  =  8 
is  consistent  with  the  observed  maximum  improvement 
in  the  value  of  it,,,  relative  to  (rr;’,  -I-  fT-j  observed 

in  Fig.  3(a).  However,  it  is  not  consistent  with  Eq.  (10) 
and  the  accompanying  discussion  of  the  preceding 
section,  which  gave  reason  for  expecting  p,„  to  increase 
monotonically  with  increasing  n. 
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Fit;.  3.  Lidar  signal  return  data  from  diffusely  reflecting  target  at 
a  2.7-km  range  using  direct  detection  during  the  period  of  essentially 
constant  average  atmospheric  extinction:  (al  standard  deviation  of 
returns  from  lasers  I  and  2  and  their  ratio  as  functions  <if  the  number 
of  pulses  averaged;  (h)  temporal  autocorrelation  and  cross  correlation 
of  laser  return  signals  as  functions  of  time  delay;  and  tc)  comparison 
of  measured  and  calculated  value  of  signal-averaged  cross-correlation 
coefficient  as  a  function  of  the  number  of  pulses  averaged. 

To  explain  this  apparent  contradiction  it  should  be 
noted  that  the  deivation  of  Eq.  (9)  assumed  Pjx  =  Pjy 
As  seen  in  Fig.  :i(b)  this  is  not  strictly  accurate;  more 
important,  it  can  be  observed  that  there  is  a  significant 
amount  of  scatter  or  noise  in  the  correlation  function 
data.  This  scatter  in  denoted  tj,  (i  =  1,2,3  =  x,y,xy, 
respectively),  appears  to  be  of  the  order  of  ±0.02  and  is 
at  least  as  large  as  p^^y  as  for  pjx  and  pjy.  When  the 
value  of  the  correlation  functions  decreases  to  values 
comparable  with  or  smaller  than  <y,,  which  occurs  of  the 
order  of  1  sec  for  the  data  shown  in  Fig.  3(b),  the  ap- 
jiarent  effect  is  to  bias  the  ratio, 

0',.. 

toward  a  smaller  value  relative  to  that  for  r  =  0.  While 
the  reason  for  this  is  not  fully  understood,  it  is  believed 
to  be  due  to  the  increa.sed  probability  of  obtaining  a 
negative  value  for  delay  times  corresponding  to  j  values 
for  which  pj,  <  r,,.  For  the  case  of  varying  atmospheric 
extinction  this  was  not  observed,  since,  as  seen  in  Fig. 
2(1)),  Pi,  >  („  over  the  entire  range  of  delay  times  con¬ 


sidered  (to  ~3  min).  Hence  the  expected  monotonic 
increase  of  p„c  with  increasing  n  was  observed  for  that 
case.  It  is  improbable  that  the  observed  maximum  is 
an  experiment  artifact  since  similar  results  have  been 
obtained  in  several  experiments  performed  at  different 
times  and  under  differing  conditions. 

3.  Joint  considerations.  Although  significant 
differences  in  the  behavior  of  the  measured  parameters 
due  to  the  presence  or  absence  of  atmospheric  drift  are 
observed  in  the  two  experiments,  there  are  noteworthy 
features  common  to  both  experiments. 

First,  agreement  between  the  measured  and  calcu¬ 
lated  values  of  pnc  vs  n  is  obtained  in  both  experiments. 
That  this  agreement  exists  despite  the  different  func¬ 
tional  dependence  on  n  in  the  two  cases  represents  a 
strong  validation  of  Eq.  (8)  and  of  the  assumptions 
made  in  its  derivation,  as  given  in  Appendix  A.  Second, 
the  relationship  among  <i„i,  a„y,  <t„(,  and  p„c  is  accu¬ 
rately  described  by  Eq.  (5)  for  both  sets  of  experiments. 
This  indicates  that  Eq.  (5)  is  valid  for  (Tnx  <^ny  <  0-2. 
The  subject  of  the  validity  of  Eq.  (5)  for  larger  values  of 
a„x  and  o^y  will  be  addressed  more  fully  below. 

Finally,  and  perhaps  most  important,  it  should  be 
noted  that  the  effect  of  using  a  dual  laser  in  combination 
with  the  signal  averaging  pr(x;ess  yielded  comparable 
values  of  in  both  cases.  For  the  case  in  which  the 
extinction  coefficient  remains  constant,  the  temporal 
correlation  of  the  individual  laser  returns  more  rapidly 
approached  zero,  the  improvement  due  to  the  p„c  term 
is  relatively  small,  and  the  decreasing  values  of  (j„(,  with 
increasing  n  are  dominated  by  the  effect  of  signal  av¬ 
eraging  on  the  individual  returns.  This  is  not  true  when 
the  atmospheric  extinction  slowly  changes  during  the 
experimental  periixl.  In  that  case,  as  seen  in  Fig.  2(a), 
the  effectiveness  of  signal  averaging  in  decreasing  the 
standard  deviation  of  the  mean  of  the  individual  laser 
return  signals  is  severely  reduced.  However,  this  is 
accompanied  by  positive  long-term  cross  correlation, 
which  results  in  the  increasing  value  of  p,„.  with  in¬ 
creasing  n  shown  in  Fig.  2(c)  and  which  effectively 
cancels  the  effect  of  atmospheric  drift  on  the  signal- 
averaged  ratio  of  the  returns,  in  accord  with  Flq.  (8).  A 
dual-la.ser  system  capable  of  obtaining  signal-averaged 
ratio  values  may,  therefore,  play  an  extremely  valuable 
role  in  a  direct-detection  DIAL  system  studying  re¬ 
flections  from  a  diffusely  reflecting  target.  Its  use  can 
result  in  a  3()-,50%  advantage  when  the  atmosphere  is 
essentially  constant  during  a  measurement.  More 
important,  it  effectively  overcomes  the  much  greater 
measurement  uncertainty  that  can  occur  in  a  slowly 
changing  atmospheric  environment. 

2.  Larfiv  Standard  Ih'cialion  and  l.arpv  Cross 
Correlation 

It  has  been  observed  that  the  standard  deviation  of 
the  backscattered  radiation  a  is  always  large  (>1 )  when 
a  heterodyne  detection  system  is  used.''’  'I'he  corre¬ 
sponding  cross-correlation  coefficient  p,  is  large  when 
the  lidar  signals  are  backscattered  from  a  specular  tar¬ 
get.  However,  for  heterodyne  detection  and  a  diffusely 
reflecting  target,  p,  is  large  (>0.9)  only  when  the  ra- 
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diation  frequencies  v  and  i>'  from  the  two  lasers  are  the 
same.  When  i'  i>',  Pc  falls  off  rapidly,  and  it  is  found 
that  p,.  <  0.05  for  frequency  separations  as  small  as  2 
cm~',  the  separation  between  adjacent  CO2  laser  lines. 
This  loss  of  cross  correlation  is  due  to  the  surface 
roughness  of  the  target,  which  gives  rise  to  different 
speckle  patterns  for  the  two  laser  beams  when  i'  ^  v'. 
The  relationship  of  cross  correlation  with  surface 
roughness  and  wavelength  is  given  by-*"  p,.  ~ 
expj— —  I'')]"!,  where  is  the  standard  de¬ 
viation  of  the  surface  roughness.  We  find  p,  <  0.05  for 
(j/  —  //)  ~  2  cm“ which  indicates  >  2  mm.  The  case 
where  v  ^  v'  will  be  discussed  in  the  following  sec¬ 
tion. 

As  an  example  of  results  for  large  standard  deviation 
(ff  ~  1.0)  and  large  cross  correlation  (p<  >  0.9),  we  will 
consider  measurements  made  of  the  backscattered  ra¬ 
diation  from  a  diffusely  reflecting  target  at  a  range  of 
2.7  km,  with  detection  in  the  heterodyne  mode  and  lx»th 
lasers  of  the  dual-laser  .system  radiating  on  the  10.6-pm 
P(20)  CO2  laser  line;  i.e.,  i'  =  Under  these  conditions, 
(T^  and  fly  were  slightly  greater  than  unity  and  p,.  ^ 
0.96. 

The  values  obtained  for  ff„v,  n,,,  ,  and  (t„c  as  functions 
of  n  for  these  experimental  conditions  are  given  in  Fig. 
4(a).  The  results  are  seen  to  differ  sharply  from  those 
given  in  Figs.  2(a)  and  9(a).  The  reduction  of  and 
ff„v  with  increasing  n  is  considerably  slower  than 
However,  tbe  high  cross-correlation  coefficient  results 
in  a  sharp  reduction  of  the  normalized  standard  devia¬ 
tion  of  the  pulse  ratio  (Tj  to  0.58.  Furthermore,  it  is  seen 
that  the  reduction  of  (t„c  with  increasing  n  is  initially 
more  rapid  than  n  when  one  properly  averages  each 
set  of  backscattered  lidar  returns  over  n  pulses  prior  to 
taking  their  ratio.  This  remarkable  result  may  result 
from  a  combination  of  three  factors:  ( I )  small  temporal 
autocorrelation  of  pulse-pair  ratio;  (2)  increa.sing 
cross-correlation  coefficient  p,„  with  increasing  n;  and 
(3)  loss  of  validity  of  Eq.  (5)  for  large  fr„,  and  <t„v. 
These  factors  are  considered  in  greater  detail  in  the 
following  subsections. 

1.  Temporal  autocorrelation  of  pulse-pair  ratio. 

The  high  cross  correlation  between  lasers  1  and  2  results 
in  a  strongly  reduced  value  of  the  temporal  autocorre¬ 
lation  of  the  ratio  of  the  pulse  pairs  (n  =  1 )  from  the  two 
lasers  P;u/,i  relative  to  those  of  the  individual  lasers, 
Pjx  and  Pjy,  for  short  delay  times.  'I'his  is  seen  in  Fig. 
4(b).  which  gives  pjx,  Pis.  and  /),(,/,)  as  functions  of 
delay  time  (}  r)  for  j  =  1, 2,4,8, . . . ,  512.  'Phis  small 
pulse-pair  autocorrelation  factor  alone  is  sufficient  to 
give  rise  to  a  standard  deviation  for  the  ratio  which 
initially  decreases  almost  as  as  seen  from  the 

curve  labeled  ratio  first,  then  averaged  in  K'ig.  4(a).  As 
noted  previously,  taking  the  ratio  of  the  individual 
pulses  prior  to  averaging  <iver  the  n  pulses  biases  the 
result.  For  the  case  of  high  cross  correlation,  the  bias 
is  relatively  small,  but.  as  seen  in  Fig.  4(a),  its  effect  is 
highly  significant. 

2.  Cross-correlat ion  effect.  When  one  properly 
averages  the  lidar  returns  (rom  the  individual  lasers  over 
the  n  pulses  prior  to  taking  their  ratio,  one  find-i  that  the 
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Fig.  4.  l.idar  signal  return  data  from  diffusely  reflecting  target  at 
a  2.7-km  range  using  heterodyne  detection  with  both  CO2  lasers  ra¬ 
diating  on  the  10.6-gm  f’(2())  laser  transition  (X  =  X'):  (a)  Standard 
deviation  of  returns  from  lasers  1  and  2  and  their  ratio  as  functions 
of  the  number  of  pulses  averaged,  showing  the  variation  for 

small  n  when  the  signals  are  properly  averaged  prior  to  taking  ratios; 
(b)  temporal  autocorrelation  of  individual  laser  return  signals  and 
of  their  ratio  as  functions  of  time  delay;  and  (c)  measured  signal- 
averaged  cross-correlation  coefficient  variation  with  number  of  pulses 
averaged  showing  the  near  approach  to  unity  with  increasing  n. 

cross-correlation  coefficient  /o„r,  as  defined  by  Eq.  (9), 
increases  with  increasing  n,  resulting  in  a  further  de¬ 
crease  of  with  increasing  n.  This  is  shown  in  Fig. 
4(c)  where  Pn,  is  seen  to  increase  from  0.964  for  n  =  1 
to  Pn,  >  0.98  for  n  >  16. 

A  proportionately  greater  increase  in  p„<.  with  in¬ 
creasing  n,  at  least  initially,  was  observed  under  the 
direct  detection  of  experimental  conditions  discussed 
previously,  with  less  dramatic  effect  on  The 

greater  effect  on  (t„j  observed  in  this  case  is  due  to  the 
closeness  of  p,„  to  unity.  This  can  be  seen  ratber  sim¬ 
ply  by  considering  the  case  where  =  a,,^  .  Then 
from  Eq.  (5)  a„^  =  “  Pm  )''’  and 

d«„i  _  _  n,iJ\  2 

dpn,  (l-p,„  **'* 

As  seen  from  F)q.  (11)  the  decrease  of  it,, j  with  increasing 
Pn,  increases  dramatically  when  approaches 
unity. 
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3.  Validity  Limit  of  Eq.  (5).  As  noted  previously,  the 
derivation  of  Eqs.  (3)  and  (5)  involved  a  first-order 
Taylor  series  expansion,  and  the  validity  of  these 
equations  is,  therefore,  limited  to  cases  where  the 
higher-order  terms  can  be  ignored.  It  is  shown  in  Ap¬ 
pendix  B  that  this  requires  a'iy  «  1. 

Since  all  the  terms  in  Eq.  (3)  have  been  measured 
independently,  it  is  possible  to  make  a  direct  experi¬ 
mental  comparison  between  the  measured  value 
(T„f  (meas)  imd  Ont  (calc)  as  calculated  on  the  basis  of  Eki. 
(5).  The  results  are  shown  in  Fig.  5  along  with  a„y.  It 
is  seen  that  for  small  n,  where  a„y  is  large,  <r„((meas)  is 
significantly  greater  than  the  value  predicted  by  Eq.  (5). 
With  increasing  n  (and  decreasing  ff„y),  ff„{(meas) 
decreases  rapidly  and  approaches  the  calculated  value. 
For  n  =  8,  where  <t„£  =  0.5,  <T„f  (meas)  «=  ff„£(calc,  and 
for  n  >  8,  the  measured  value  of  (r„(  remains  in  close 
agreement  with  the  value  predicted  by  Eq.  (6).  This 
result  indicates  that  Eq.  (5)  represents  a  good  approx¬ 
imation  for  <  0.25. 

In  general,  we  have  always  observed  <r„£(meas)  » 
fT„j(calc)  when  (T„y  »  0.5.  In  the  signal  averaging 
process,  this  large  a  value  of  has  generally  been 
limited  to  small  values  of  n  in  our  experiments  (i.e.,  to 
n  <  8  for  the  example  in  Fig.  4(a)j.  Therefore,  in  this 
region  of  small  n,  there  may  be  a  relatively  rapid  de¬ 
crease  in  the  measured  value  of  with  increasing  n 
as  the  value  goes  from  (T„((meas)  »  <7„{(calc)  for  n  = 
1  to  (Tnjlmeas)  =«=  (T„£(calc)  for  n  =  8.  Indeed,  in  Figs. 
4(a)  and  (5),  the  region  where  <T„{(meas)  decreases  with 
increasing  n  more  rapidly  than  jg  limited  to  n  < 
8.  Furthermore,  it  should  be  noted  that  (calc)  varies 
everywhere  more  slowly  than 

Large  Standard  Deviation  and  Small  Cross 
Correlation 

As  noted  above,  the  cross  correlation  is  extremely 
small  for  heterodyne  detection  and  backscatter  from  a 
diffusely  reflecting  target  when  the  two  lasers  are  op¬ 
erating  at  different  frequencies.  For  frequency  dif¬ 
ferences,  I’  —  v',  as  small  as  2  cm“',  pc  is  typically  near 
zero  ( <0.05).  The  value  of  p„,  tends  to  increase  slowly 
with  increasing  n,  but  even  after  averaging  over  512 
pulses  Pne  values  of  ~0.25  ±0.1  are  typical.  Since  dif¬ 
ferential  absorption  experiments  require  v  v' ,  this 
small  cross  correlation  represents  a  serious  limitation 
to  the  effectiveness  of  using  lidar  backscatter  ratios  for 
reducing  the  effective  standard  deviation  of  heterodyne 
DIAL  measurements. 

In  addition,  the  combination  of  large  fluctuations  and 
near-zero  cro.ss  correlation  between  the  lidar  returns  at 
the  two  frequencies  leads  to  a  strong  bias  in  their  ratio 
and  a  large  excess  variance.  As  a  result,  the  measured 
standard  deviation  of  the  averaged  ratio  initially  (for 
small  n )  is  greater  than  would  be  obtained  by  simply 
considering  the  standard  deviation  of  the  individueil  sets 
of  lidar  returns  with  zero  cross  correlation,  i.e.,  (t„  t> 

This  is  shown  in  F'ig.  6,  which  gives  the  values  t)f 
fTnjImeas),  (T^jlcalc),  rTni,and  fT„y, alias 

functions  of  n,  obtained  experimentally  with  hetero- 
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Fig.  5.  Comparison  of  measured  and  calculated  values  of  a„(  as  a 
function  of  the  number  of  pulses  averaged  based  on  the  data  used  in 
Fig.  4.  Agreement  is  seen  to  be  excellent  for  a„y  <  0.5. 
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Fig.  6.  Standard  deviation  as  functions  of  the  number  of  pulses 
averaged  for  lidar  signal  returns  from  the  diffusely  reflecting  target 
at  a  2.7-km  range  using  heterodyne  detection  with  lasers  I  and  2  ra¬ 
diating  on  the  P(20)  and  P(22)  COj  laser  transitions  of  the  10.6-pm 
hand,  respectively  (A  X'). 


dyne  detection  of  backscattered  lidar  signals  from  a 
diffusely  reflecting  target  with  e  y  v'.  Specifically,  for 
the  results  shown  in  Fig.  6,  a  flame-sprayed  aluminum 
plate  at  a  2.7-km  range  was  used  as  the  diffusely  re¬ 
flecting  target,  with  lasers  1  and  2  operating  on  the 
10.6-pm  P(20)  and  P(22)  CO2  laser  transitions,  re¬ 
spectively.  For  these  lines,  c  -  P  =  1.8  cm"'.  From  the 
figure,  it  can  be  seen  that  use  of  backscattered  signal 
ratios  has  some  limited  value  in  reducing  measurement 
uncertainty  at  the  larger  n  values  considered  but  exacts 
a  significant  penalty  at  small  n  because  of  the  bias  ef¬ 
fect.  It  should  be  noted  that  the  results  shown  are 
consistent  with  the  conclusions  reached  in  the  preceding 
section  regarding  the  limits  of  the  validity  of  Eq.  (5). 

V.  Conclusions 

The  previously  derived  relationships  expressed  by 
Eqs.  (3)  and  (4)  have  been  extended  to  encompa.ss  the 
effects  of  both  signal  averaging  and  cross  correlation  on 
the  variance  of  the  ratio  of  averaged  backscattered 
signals  from  the  two  lasers  of  a  dual-laser  DIAL  system. 
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This  extension  involved  deriving  the  relationship  of  the 
signal-averaged  cross-correlation  coefficient  with  the 
temporal  autocorrelation  and  cross  correlation  of  the 
individual  pulses. 

Extensive  experimental  results  using  both  heterodyne 
and  direct-detection  CO)  lidar  measurements  were 
presented  and  shown  to  be  consistent  with  the  derived 
relationships.  The  DIAL  experiments  investigated  the 
effectiveness  of  using  the  ratio  of  signal  returns  from  a 
dual-laser  system  in  conjunction  with  signal  averaging 
in  improving  measurement  accuracy.  In  general,  sig¬ 
nificant  improvement  was  obtained  through  such  a 
combination  of  signal  averaging  and  then  taking  ratios. 
However,  the  degree  of  improvement  was  strongly  de¬ 
pendent  on  the  relative  values  of  the  standard  deviation 
of  the  lidar  signals  and  the  temporal  cross  correla¬ 
tion. 

Specifically,  with  diffuse  targets  and  direct  detection 
(fT„j  and  Pn,  small),  improvement  due  to  cross  correla¬ 
tion  was  slight  when  the  atmospheric  extinction  was 
essentially  constant  throughout  the  measurement  pe¬ 
riod;  therefore,  under  these  conditions,  use  of  a  dual¬ 
laser  system  does  not  offer  significant  improvement  in 
DIAL  measurement  accuracy.  However,  when  the  at¬ 
mospheric  extinction  is  slowly  but  constantly  varying 
during  the  measurement  period,  use  of  ratios  with  a 
dual-laser  system  effectively  eliminates  this  variation 
as  a  major  source  of  error. 

For  heterodyne  detection  of  returns  from  a  diffusely 
reflecting  target  it  was  found  that  both  the  cross-cor¬ 
relation  and  the  standard  deviation  are  initially  large 
when,  and  only  when,  c  =  i''.  On  signal  averaging,  ffnj 
decreases  with  increasing  n  while  Pn,  remains  large  or 
increases.  Under  these  conditions,  the  combination  of 
signal  averaging  and  then  taking  ratios  was  shown  to  be 
most  effectise  with  an  initial  decrease  of  more  rapid 
than  n  ~ '  -  possible.  On  the  other  hand,  for  heterodyne 
detection  with  backseat  ter  from  a  diffusely  reflecting 
target  and  v  r'.  p,  is  small,  and  the  standard  deviation 
of  the  individual  laser  returns  is  large.  This  is  the 
condition  for  maximum  biasing  of  the  ratio  and  results 
in  a  measured  value  of  which  is  greater  than  would 
he  obtained  by  considering  the  lidar  signal  returns  of  the 
individual  lasers.  'I'herefore,  use  of  heteodyne  detec¬ 
tion  in  a  differential  absorption  measurement  requires 
signal  averaging  over  a  large  number  of  pulses  (n  >  10) 
at  each  wavelength  prior  to  taking  the  ratio  to  avoid 
having  this  bias  effect  introduce  excess  error.  However, 
even  after  signal  averaging,  the  effective  cross  correla¬ 
tion  p,„  remains  small,  and,  therefore,  use  of  a  dual-laser 
system  was  luund  to  be  of  limited  value  under  these 
experimental  conditions. 

In  addition,  the  limits  of  validity  of  Kqs.  (O)  and  (5) 
were  investigated  by  carrying  out  a  'I'aylor  series  ex¬ 
pansion  of  the  ratio  to  the  second  order.  It  was  found 
that  higher-order  terms,  which  were  ignored  in  deriving 
the  relationship,  could  make  a  significant  contribution 
for  large  values  of  .  'This  agreed  with  expi  i  imentally 
observed  differences  between  <T,j,(meas)  and  rrUfcalc) 
lor  IT,,  >  O.'J,").  where  fi;’,(calcl  was  determined  on  the 
basis  of  Fk|.  (’>).  However,  for  it,,,  <  0.2*).  generally  g(M)d 
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agreement  between  the  experimental  and  calculated 
values  of  a^^was  obtained. 

The  Lincoln  Laboratory  portion  of  this  work  was 
supported  by  the  National  Aeronautical  and  Space 
Administration  and  the  Air  Force  Engineering  and 
Services  Center.  The  University  of  Maryland  portion 
of  this  work  was  supported  by  the  Department  of  En¬ 
ergy. 

Appendix  A 

To  derive  Eqs.  (5)-(7),  which  describe  the  effect  of 
signal  averaging  on  cross  correlation  and  on  the  variance 
of  the  lidar  return  ratios,  several  quantities  must  be 
defined.  We  first  define  Ikx  =  Ix(tk)  and  hy  =  + 

A<)  as  the  normalized  deviation  of  the  feth  lidar  pulse 
return  signals,  Pkx  and  P*y,  occurring  at  tinrms  and 
4  +  A<,  respectively,  from  their  mean  values  Px  and  Py  \ 
i.e., 

/*)  =  -  T\  )/?.  and  Uy  =  (H^y  -  Py)/Vy,  (At) 

where,  as  in  Eq.  (,‘l),  subscripts  x  and  y  refer  to  lasers  1 
and  2,  respectively.  The  normalized  variance  of  the  full 
set  of  pulses  from  laser  1  is  then  defined  as 

0:  =  =:^  L  /L.  (A2| 

1  *=i 

where  F  is  the  total  number  of  pulses  from  each  laser  in 
the  data  set.  (F  =  6144  in  most  experiments  described 
in  this  paper.)  The  corresponding  temporal  autocor¬ 
relation  coefficient  is  given  by 

f>ii  =  +  Jt\) 

n\ 

1  y -I 

~  ..  L  (A,tl 

fT;(l  -y)  *  =  1 

where  t  is  the  time  interval  between  pulses  from  the 
individual  laser.  The  definitions  of  the  variance  and 
temporal  autocorrelation  coefficients  for  laser  2,  nf  and 
P;y ,  are  essentially  the  same  as  tho.se  given  in  Eqs.  (A2) 
and  (AD)  with  X  ♦  y.  The  pulse-pair  cross-correlation 
coefficient  p,  is  defined  as  the  covariance  of  the  paired 
lidar  returns  from  lasers  1  and  2  normalized  to  their 
standard  deviation;  that  is, 

_  IT,'  _  -I-  An> 

fT ,  (7 ,  fT ,  rr  ^ 

=  — L  /*,/»'  (A4) 

IT, IT,  I 

When  considering  temporal  cross  correlation  for  delay 
times  of  the  order  of  jr  (i.e.,;)-  ±  A)  where  At  «  jr). 
one  has 

_  ■  /, \u  I/,  III.  +  At  +  it)) 


= - :: - 1/h./.i..;k  iA-^i 

n  ,  11 '  I  I  -  1 1 

Using  the  above  definitions,  one  may  evaluate  the 
variance  and  the  cross  correlation  of  the  F  lidar  returns 
from  both  lasers  after  they  have  been  averaged  over  n 
returns.  'Fhis  calculation  has  already  been  carried  out 
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for  the  variance  of  the  individual  lasers  in  Ref.  19,  where 
it  was  shown  that 

f  n-I 

=  (af/n)  1  +  2  L  d  “  JMpj,  ,  (A6) 

I  y-i 

where  i  = 

The  corresponding  averaged  cross-correlation  coef¬ 
ficient  is  given  by 


p-«- 
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■  +  fny\ 

n 

1 
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\  n  A  n  I 


//|r-(n-ll|,  + 

.  +  /r. 

i 

(A7) 
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(A8) 

The  terms  within  the  first  bold  parentheses  of  Eq. 
(A8)  can  be  broken  int<»  pairs  which  refer  to  the  same 
subset  of  n  pulses;  i.e.,  k  =  I  to  n,  n  -+  1  to  2n, . . .  I"  — 
(n  —  1)  to  r,  as  shown.  The  time  delay  differs  in  the 
two  terms  of  such  a  pair  in  that  the  delay  will  be  (j  t  + 
Af )  for  one  term  and  ijr  -  At)  for  the  other  term.  In 
our  experiments,  ;  t  is  some  multiple  of  100  msec,  while 
At  ^  5-10  ^sec.  The  atmosphere  can  be  considered 
frozen  over  the  time  delay  difference  2At ,  and  this  dif¬ 
ference  can,  therefore,  be  ignored.  In  addition,  since 
both  Ikx  and  /*,  are  normalized,  one  may,  as  a  reason¬ 
able  approximation,  assume  that  when  averaged  over 
the  full  set,  Ikxhk+j)\  ^  hk+j)xlky  In  that  ca.se,  the  two 
terms  dealing  with  the  same  pulse  subset  may  be  con¬ 
sidered  equivalent,  and  from  Eqs.  (A4)  and  (A5) 
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nl  (7  ttx 


.n  —  1\ 

I  +  2U  -  -  /Miv  + 


+  2ir-(n  -  I)|t7^fT, 


or 


/»u.  = 


f>.  +  L  ( 1 

/“I 


-  ;/!')( I  -  )/n)fij„ 


(A9I 


(AI0» 


Equations  ( A91  and  ( AlO)  assume  equal  weighting  of 
interpolated  missing  terms  in  Eq.  (A8),  as  discussed  in 
Ref.  19.  The  derivation  of  Eq.  (AlO)  closely  parallels 
the  derivation  of  <i^  given  there  for  the  segmentally 
averaged  ca.se.  It  was  shown  that  for  a  sufficiently  large 
sample,  the  (1  -y/I'lterm  »  1,  and  Eq.  (AlO)  effectively 
reduces  to  the  form 


fT,  ff, 

n  fJfj  I 


P,  +  2  (1  -  j/n>p,„ 


(All) 


In  a  similar  fashion,  one  may  use  the  approach  de¬ 
scribed  in  Appendix  C  of  Ref.  19  to  show  that  Eq.  ( A1 1) 


is  obtained  directly  when  a  running  average  rather  than 
a  segmental  average  is  used  to  define  the  mean  value  of 
the  n  pulses  being  integrated. 

Appendix  B:  Calculation  of  Variance  of  ^(=x/y)  to 
Second  Order 

In  general,  the  Taylor  expansion  of  a  function  z  of  two 
variables,  x  and  y ,  evaluated  at  a  point  p  corresponding 
to  the  averaged  values,  x  and  y,  is  given  to  second  order 
in  the  derivatives  by 

z=l(x.y)=f(x,y) 


+  ^l 

dx 


t>f 

(l-i-)  +  —  (y-v) 

p  dy  p 


dy\ 


(y-y)2+2^  (l-l)(j--y) 

P  dxay  p 


(Bl) 


while  the  average  value  z  is 
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(B2) 


The  variance  of  z  may  be  given  by 


a]  =  ({z  -  z)“> 
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(B.l) 


where  all  the  derivatives  are  evaluated  at  point  p  =  x,y. 
For  z  =  ^  =  x/y 
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Substituting  into  Eq.  (B3)  and  normalizing, 
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(x2/y2) 


(X  -  x)'^  2(x  -  xMv  -  vl 


(v  -  v)^  f(j  -  X  >-'(>  -  V) 

r  i  x-v 

2(x  -  x)(,v  -  \ )-'  ^  (  V  -  yi  'j 


f(x  -  ?»-'(v  -  V)-  2(x  -  x)(v  -  v)'* 
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(B4) 


In  Eq.  (B4),  rr,  and  tTy  represent  the  normalized  stan¬ 
dard  deviation  values.  The  equation  may  be  re¬ 
written 


—  <Ty((T^  +  <t:;  —  ) 


2(x  -  xl(  v  -  y)| 
.XV  I 


(B-S) 


The  term  in  braces  on  the  rljjht-hand  side  of  Eq.  (B5) 
is  just  the  solution  obtained  for  the  variance  of  the  lidar 
return  of  tr'^  for  the  first-order  Taylor  series  expansion 
and  as  given  in  Eq.  (2).  In  that  equation,  tr'i  is  identified 
with 


((V  -  ,vl-> 


It  is  apparent  that,  for  values  of  a;  (or  in  the  signal 
averaged  case)  approaching  unity,  higher-order  terms 
in  the  Taylor  series  expansion  play  a  significant  role  and 
cannot  be  ignored.  Therefore,  the  validity  of  Eqs.  (3) 
and  (5)  is  limited  to  the  cases  where  fr'i  «  1  or  «  1, 
respectively.  This  is  consistent  with  the  experimental 
results  given  in  the  paper,  which  indicate  serious  dis¬ 
crepancies  between  measured  values  of  (t^j,  and  values 
calculated  on  tbe  basis  of  Eqs.  (3)  or  (5)  occur  when  c'iy 
>  0.25. 

Although  the  moments  of  each  of  the  third-  and 
fourth-order  terms  in  Eq.  (B4)  can  be  calculated  when 
the  distribution  function  of  tbe  lidar  signals  is  known, 
sucb  calculations  are  beyond  tbe  scope  of  this  paper. 
However,  it  should  he  noted  that  while  the  third-order 
terms  would  be  zero  for  a  symmetrical  distribution 
function,  lidar  return  distributions  tend  to  be  .strongly 
asymmetric  (i.e.,  lognormal,  negative  exponential,  or  k 
distribution)  and  dependent  on  experimental  condi¬ 
tions.  Tberefore.  tbe  tbird-order  term  must  be  re¬ 
tained. 

It  was  observed  that  when  the  Taylor  series  expansion 
is  limited  to  first  order,  ai  =  This  is  not  true  for 

higher-order  expansions. 
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Appendix  B 

The  followiny  is  a  preprint  of  a  journal  article  to  be  published  in 
Applied  Optics  entitled  "Laser  Remote  Sensing  of  Atmospheric  Ammonia  Using 
CO2  LIDAR  System." 
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LASER  REMOTE  SENSING  OF  ATMOSPHERIC  AMMONIA  USING  A  CO2  LIDAR  SYSTEM* 


A.  P.  Forcet,  D.  K.  Kil linger,  W.  E.  DeFeo  and  N.  Menyuk 
Lincoln  Laboratory,  Massachusetts  Institute  of  Technology 
Lexington,  Massachusetts  02173-0073 


Abstract 


A  CO2  differential -absorption  LIDAR  system  has  been  used  for  the  remote 
sensing  of  ammonia  in  the  atmosphere.  For  CO2  LIDAR  returns  backscattered 
from  topographic  targets  at  ranges  up  to  2.7  km,  the  path-averaged  sensitivity 
of  the  DIAL  system  was  5  parts  per  billion  of  NH3.  Concentrations  of 
atmospheric  ammonia  were  found  to  vary  during  the  day  from  undetectable  levels 
(<  5  ppb)  to  as  high  as  20  ppb,  depending  upon  temperature  and  humidity 
conditions. 


*This  work  was  supported  by  The  National  Aeronautical  and  Space  Administration 
and  the  Air  Force  Engineering  and  Services  Center. 

tVisiting  Scientist  from  1).  S.  Army  Chemical  Research  and  Development  Center, 
Aberdeen  Proving  Ground. 
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LASER  REMOTE  SENSING  OF  ATMOSPHERIC  AMMONIA  USING  A  CO2  LIDAR  SYSTEM* 


A.  P.  Forcet,  D.  K.  Killinger,  W.  E.  OeFeo  and  N.  Menyuk 
Lincoln  Laboratory,  Massachusetts  Institute  of  Technology 
Lexington,  Massachusetts  02173-0073 


I.  Introduction 

Single-ended  di fferential -absorption  LIDAR  (DIAL)  has  been  shown  to  be  a 
sensitive  method  for  the  long-range  reaote  sensing  of  trace  molecular  species  in 
the  atmosphere.  DIAL  systems  operating  in  the  infrared  have  been  used  to  moni¬ 
tor  several  trace  atmospheric  constituents^"^  including  C2H4,  NO,  SO2, 
and  O3  and  work  is  underway  to  increase  the  number  of  atmospheric  species 
which  can  be  detected. 

In  this  paper,  the  use  of  a  pulsed  infrared  C02"l3ser  DIAL  system  for  the 
remote  sensing  of  atmospheric  NH3  is  reported.  For  CO2  LIDAR  returns 
backscattered  from  targets  at  a  range  of  2.7  km,  the  experiments  demonstrate  the 
ability  of  the  DIAL  system  to  remotely  sense  the  presence  of  NH3  with  a  path- 
averaged  detection  sensitivity  of  5  parts  per  billion  (ppb).  The  concentration 
of  ambient  atmospheric  NH3  was  observed  to  vary  between  undetectable  levels 
(<  5  ppb)  to  values  as  high  as  20  ppb.  In  addition,  an  apparent  negative 
correlation  between  the  concentration  of  atmospheric  NH3  and  relative  humidity 
or  water  vapor  aerosol  concentration  was  also  observed. 

Atmospheric  ammonia  is  produced  by  a  combination  of  natural  and 
anthropogenic  sources^.  Its  detection  and  measurement  is  important  in  that 
its  concentration  is  a  critical  parameter  in  the  rates  of  formation,  transport, 

*This  work  was  supported  by  The  National  Aeronautical  and  Space  Administration 
and  the  Air  Force  Engineering  and  Services  Center. 

tVisiting  Scientist  from  U.  S.  Army  Chemical  Research  and  Development  Center, 
Aberdeen  Proving  Ground. 
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transformation  and  removal  of  the  sulfur  and  nitrogen  aerosols  which  play  an 
important  role  in  the  generation  of  acid  rain^.  Previous  measurements  of 
atmospheric  NH3  have  generally  relied  on  wet -chemical  techniques,  which 
require  the  collection  of  atmospheric  samples^*®.  The  sensitivity  of  such 
techniques  along  with  the  period  of  ti«  required  for  sample  collection  have 
limited  the  time  resolution  of  the  measurements  to  periods  of  several  hours.  In 
contrast,  DIAL  systems  are  capable  of  in  situ  real-time  measurements,  and  may 
thus  play  a  significant  role  in  future  atmospheric  measurements.  This  paper 
presents  preliminary  experiments  which  indicate  the  capability  of  the  CO2  DIAL 
technique  for  measurement  of  atmospheric  NH3  concentrations. 

II.  Experimental  Apparatus 

A  schematic  of  the  dual-C02  laser  DIAL  system  used  in  our  experiments  is 
shown  in  Fig.  1  and  has  been  described  previously. The  pulsed  transmitted 
radiation  was  provided  by  two  line-tunable  mini -TEA  CO2  lasers.  The  CO2 
lasers,  typically,  have  a  linewidth  on  the  order  of  0.1  cm"^,  an  output  energy 
of  10  mJ/pulse,  a  pulselength  of  100  nsec,  and  operate  at  a  PRF  of  15  Hz. 

As  seen  in  Fig.  1,  the  outputs  of  the  two  lasers  were  collimated  and  then 
joined  with  a  50/50  beam  splitter.  Portions  of  each  beam  were  directed  to  a 
spectrometer,  a  pyroelectric  detector  to  normalize  the  pulse  intensities,  and 
through  a  Pyrex  absorption  cell  (33.5  cm  long  x  ?.54  cm  diameter  with  BaF2 
windows)  to  calibrate  the  NH3  absorption  levels.  The  remainder  of  the  laser 
radiation  was  sent  through  a  lOX  beam  expander  and  directed  to  either  a  topo¬ 
graphic  target  (the  wall  of  a  building),  a  f lame-sprayed  aluminum  plate  ,  or  a 
one-inch  retroref lector,  all  at  a  range  of  2.1  km.  Backscattered  LIDAR  returns 
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re  collected  with  a  30  cm  Cassegrain  telescope  and  focused  onto  a  cooled 
Idle  detector  (D*  =  2.8  x  10^0  cm  Hz^/^  W~l)*  The  signals  from  the 
tectors  were  recorded  and  analyzed  by  a  computer  based  data  acquisition  system 
ich  calculated  the  normalized  differential  absorption  of  the  laser  radiation 
'■ough  the  laboratory  absorption  cell  and  that  of  the  LIDAR  returns.  For 
libration  purposes,  a  large,  chemically-inert,  polypropylene  tank  (104  cm  long 
52.5  cm  diameter  with  15°  slanted  polyethylene  windows)  could  be  placed 
:side  the  window  in  the  path  of  the  laser  beams  to  measure  the  absorption  of 
?  laser  radiation  for  known  amounts  of  NH3  in  the  tank;  a  muffin  fan  was 
)ced  inside  the  tank  to  reduce  stratification  of  NH3  in  the  tank. 

Normally,  for  dual -laser  OIAL  operation  Laser  1  was  tuned  to  a  transition 
?quency  which  had  low  NH3  absorption  and  Laser  2  was  tuned  to  a  frequency 
ich  had  high  NH3  absorption;  the  lasers  were  separately  triggered  with  a 
ie  delay  of  SO-us  in  order  to  freeze  the  effect  of  atmosphere  turbulence 
;ween  the  two  output  pulses. ^  Under  some  experimental  conditions,  however, 
ly  one  laser  was  used  for  the  DIAL  experiments  and  was  sequentially  tuned  on 
I  of f-resonance;  this  single-laser  OIAL  measurement  was  necessary  because  one 
the  lasers  in  Fig.  1  would  not  lase  on  the  required  transitions  for  use  with 
!  strongest  NH3  absorption  line.  The  distinction  between  single  and  dual- 
,er  DIAL  measurements  will  be  made  clear  in  the  discussion  of  experimental 
a. 

.  Choice  of  Laser  Frequencies 

The  choice  of  CO^  laser  frequencies  for  the  remote  sensing  of  NH3  is 
ermined  by  several  factors.  First,  the  choice  of  frequencies  is  limited  to 
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Table  I 


CO?  Laser  Line  Wavelength  1^3  Absorption^ 

Coefficient 
A(ym)  o  (atm*cm)“^ 


*R(30) 

9.217 

56. 

*R(26) 

9.237 

0.1 

R(16) 

9.291 

12.7 

R(10) 

10.316 

0. 5 

R(8) 

10.331 

20. 

R(6) 

10.346 

26. 

*P(30) 

10.693 

0.9 

*P(32) 

10.716 

13.7 

P(34) 

10.738 

14.0 

Atmospheric 

Attenuation 


a  km"^ 


B 


Summer 

Winter 

.23 

.06 

.25 

.07 

.34 

.10 

.28 

.07 

NA 

NA 

NA 

NA 

.28 

.05 

.27 

,05 

.26 

.05 

From  References  12  and  13. 

From  Reference  14. 

On-Resonance  and  Off -Resonance  CO?  Laser 
Line  pairs  chosen  for  DIAL  measurements 
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FIGURE  CAPTIONS 


Fig.  1,  Schematic  of  the  dual -laser  LIOAR  system  used  for  the  remote  sensing 
of  ammonia. 

Fig.  2.  Time  variation  of  LIDAR  returns  of  10.69-pm  P(30)  and  10.71-wn  P(32) 
radiation  passing  through  the  tank  and  reflected  from  a  topographic 
target  at  a  range  of  2.7  km  after  injecting  0.2  cm^  of  28%  aqueous 
NH3. 

Fig.  3.  Dual  laser  DIAL  measurements  of  ambient  atmospheric  NH3  using  the 
10.69-pm  P(30)  and  10.71-pm  P(32)  CO2  laser  lines  over  a  2.7  km 
range  with  26%  relative  humidity. 

Fig.  4.  Single  laser  DIAL  measurements  of  ambient  atmospheric  NH3  using  the 
9.24-pm  R(26)  and  9.22-pm  R(30)  CO2  laser  lines  over  a  2.7  km  range 
with  26%  relative  humidity. 

Fig.  5.  DIAL  measurements  of  ambient  atmospheric  NH3  using  the  9.24-pm  R(26) 
and  9.22-pm  R(30)  CO2  laser  lines  over  a  2.7  km  range  showing  the 
effect  of  a  change  in  relative  humidity  from  54%  to  41%. 
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Ambient  concentrations  varying  between  zero  and  20  ppb  were  measured  in  real 
time,  and  the  system  was  shown  to  be  capable  of  detecting  atmsopheric  NH3 
concentration  variations  of  a  few  parts  per  billion,  in  reasonable  agreement 
with  the  calculated  sensitivity. 

The  concentration  of  ambient  NH3  was  found  to  be  negatively  correlated 
with  the  relative  humidity.  This  was  shown  to  be  consistent  with  the  changes 
that  are  known  to  take  place  in  atmospheric  aerosols  with  changes  in  relative 
humidity. At  humidities  above  SOX  the  NH3  is  washed  from  the  air  by 
the  wet  particles,  but  as  the  relative  humidity  drops  below  about  50%  NH3  is 
released  by  the  aerosols.  While  the  experimental  data  is  consistent  with  this 
hypothesis,  further  studies  will  be  required  in  order  to  further  document  the 
effect. 
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NH3  is  extremely  soluble  in  water  (251.6  q/l)  and  is  rapidly  removed  from 
the  atmosphere  by  aqueous  aerosols,  which  provides  a  reasonable  explanation  for 
the  observed  NH3  change  in  Fig.  5.  .  Aerosol  particle  growth  is  a  function 
not  of  the  quantity  of  precipitable  water  in  the  air  but  of  the  relative 
humidity,  being  slight  for  relative  humidities  less  than  50%  and  increasing 
rapidly  with  increasing  relative  humidity. This  aerosol  growth  factor  may 
also  be  correlated  to  the  fact  that  the  solubility  of  gaseous  ions  onto  the 
surface  of  atmospheric  aerosols  has  been  shown  to  depend  significantly  upon  a 
Critical  Relative  Humidity  (CRH),  being  very  low  for  values  less  than  the  CRH, 
and  extremely  soluble  for  large  values. 20  Small  changes  in  the  relative 
humidity  around  this  critical  value  can  result  in  a  very  large  change  in  the 
concentration  of  the  gaseous  species.  Typical  CRH  values  range  in  the  order  of 
50  to  70%  for  several  ionic  species  in  aqueous  solution  including  NH4NO3  and 
NH4CI.  While  the  CRH  value  for  NH3  (i.e.,  NH4OH  in  solution)  is  not  di¬ 
rectly  known,  it  is  anticipated  that  a  similar  value  would  apply.  If  this  is 
the  case,  a  reasonable  explanation  for  the  data  observed  in  Fig.  5  is  that  for  a 
relative  humidity  of  54%,  the  aerosol  particles  are  essentially  wet  and  have 
absorbed  most  of  the  soluble  NH3  gas  from  the  air;  when  the  relative  humidity 
drops  to  41%,  the  aerosol  particles  act  as  dry  particles  releasing  the  absorbed 
NH3.  Obviously,  much  further  expermental  work  will  be  required  in  order  to 
determine  the  validity  of  this  assertion  as  well  as  the  underlying  rate 
constants. 

VI.  Conclusions 

A  CO2  laser  DIAL  system  has  been  used  to  remotely  detect  the 
path-averaged  concentration  of  NH3  in  the  atmosphere  over  a  range  of  2.7  km. 


produced  data  that  varied  by  only  about  10%  in  the  determination  of  Pa/Pa '  in 
Eq.  (1);  this  factor  corresponds  to  an  uncertainly  in  the  determination  of  the 
NH3  concentration  of  about  3  ppb. 

The  resultant  DIAL  data  are  given  fn  Fig.  4  and  show  a  path-averaged  NH3 
concentration  of  15  ±  5  ppb  over  a  period  of  about  one  hour;  these  measurements 
were  made  when  the  relative  humidity  was  about  26%.  As  seen  in  the  data,  the 
overall  accuracy  of  the  DIAL  measurement  is  better  than  that  shown  in  Fig.  3  due 
to  the  use  of  the  stronger  absorption  line,  R(30). 

The  single-laser  DIAL  experiments  were  repeated  on  another  day;  in  this 
case,  however,  the  relative  humidity  was  initially  54%  (temperature 
approximately  50  OF)  and  there  was  a  considerable  amount  of  haze.  The  data 
are  given  in  Fig.  5  and  indicate  essentially  zero  NH3  concentration  present  in 
the  atmosphere  during  the  morning  hours  when  the  humidity  was  high.  Early  in 
the  afternoon  the  haze  burned  off  and  there  was  a  drop  in  humidity  to  41% 
(temperature  approximately  60°F).  With  this  change  in  humidity,  NH3 
concentration  rose  to  20  i  5  ppb.  This  negative  correlation  between  relative 
humidity  and  measured  NH3  concentration  was  also  observed  on  several  other 
days  with  similar  meteorological  behavior.  This  observed  NH3  and  humidity 
anti-correlation  is  not  believed  to  be  caused  by  changes  in  the  concentration  of 
water  vapor  leading  to  a  water  vapor  induced  differential-absorption  signal.  In 
this  regard,  it  is  estimated  that  the  absolute  water  vapor  concentration  varied 
by  8%  from  5.0  Torr  (54%  relative  humidity,  50OF)  to  5.4  Torr  (41%  relative 
humidity,  6OOF).  Using  this  variation  in  the  (a  -  a')  term  in  Eq.  (1),  one 
estimates  that  such  a  change  would  only  influence  the  deduced  concentration  of 
NH3  by  0.3  ppb,  a  value  negligible  in  comparison  to  the  observed  NH3  change. 
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very  low  output  energy  {<  0.5  mJ/pulse)  of  the  CO2  laser  on  the  R(30)  or  R(26) 
line.  Details  of  these  two  experiments  follow. 

In  the  first  experiment  Laser  1  was  tuned  to  the  P(30)  reference  line  and 
Laser  2  was  tuned  to  the  P(32)  line  near  9.7  um.  The  LIDAR  returns  from  the  two 
lasers  were  collected  and  averaged  for  500  pulses  (30  seconds).  The  deduced 
path-averaged  NH3  concentration  Is  plotted  against  the  time  of  day  In  Fig.  3. 
Peaks  of  10  to  20  ppb  were  occasionally  observed  above  a  baseline  scatter  of 
10  ppb;  this  Includes  an  estimated  bias  of  8  ppb  due  to  the  (o  -  o')  term  In 
Eq.  (1).  While  the  observed  peaks  are  not  much  higher  than  the  estimated  back¬ 
ground,  the  data  Is  suggestive  that  the  observed  DIAL  signal  is  due  to 
atmospheric  NH3,  and  not  an  artifact,  since  other  atmospheric  species  which 
absorb  radiation  at  the  P(32)  laser  line,  such  as  CHCI3  and  CF2CI2,  also 
have  similar  absorbance  at  the  P(30)  llne.^^  In  addition,  the  ID  to  20  ppb  of 
NH3  observed  Is  In  good  agreement  with  the  2  to  20  ppb  that  Is  considered  to 
be  normal  in  the  atmosphere.^® 

Measurements  using  the  more  optimal  R(30)  and  R(26)  lines  near  9.2ym  band 
were  attempted  to  take  advantage  of  greater  sensitivity  available  at  these 
frequencies.  Because  the  CO2  laser  energy  was  limited  to  about  0.5  mJ  on  the 
R(30)  or  R{26)  line,  the  returns  from  a  diffuse  target  at  a  range  of  2.7  km 
produced  a  signal-to-nolse  ratio  of  less  than  one.  To  increase  the 
signal-to-nolse  ratio,  the  LIDAR  was  directed  to  a  one-inch  ret roref lector  at 
the  same  range.  Although  this  gave  a  stronger  return,  the  small  target  size  and 
the  deflection  of  the  beam  by  atmospheric  turbulence  increased  the  standard 
deviation  of  the  individual  returns.  However,  signal  averaging  over  500  pulses 


Initial  DIAL  calibration  measureiients  of  NH3  were  made  using  the  poly¬ 
propylene  tank  shown  in  Fig.  1.  The  tank  was  purged  with  dry  N2  for  one  hour 
and  placed  in  the  LIDAR  path.  The  100%  transmittance  baselines  for  the  P(30) 
and  P(32)  lines  were  obtained  by  averaging  the  normalized  LIDAR  returns  from  a 
building  at  a  range  of  2.7  km  for  1000  pulses.  A  0.2  cm3  sample  of  28% 
aqueous  NH3  was  then  injected  into  the  tank.  The  change  in  the  LIDAR  returns 
as  the  NH3  solution  evaporated  in  the  tank  is  shown  in  Fig.  2. 

As  seen  in  Fig.  2,  the  transmittance  values  of  the  P(30)  and  P(32)  lines, 
after  a  steady  state  was  reached,  were  0.81  and  0.12,  respectively.  The  concen¬ 
tration  of  NH3  in  the  tank  was  estimated  to  be  220  ppm;  the  0.2  cm3  sample 
of  28%  NH3  with  a  specific  gravity  of  0.898  g/cm3  yields  1.8  x  10^^  mole¬ 
cules  of  NH3  in  the  319  liter  tank.  With  these  transmission  values,  NH3 
concentration  levels,  and  a  round-trip  path  length  through  the  cell  of  208  cm, 
the  experimental  absorption  coefficients  for  the  P(30)  and  P(32)  lines  were  1.1 
(atm  cm)"^  and  12  (atm  cm)"^,  respectively.  These  values  are  in  good 
agreement  with  those  measured  in  the  laboratory  absorption  cell  and  the  data 
given  in  Table  I. 

V.  Laser  Remote  Sensing  of  Atmospheric  NH3 

The  DIAL  system  was  used  to  detect  and  measure  the  path-averaged 
concentration  of  ambient  atmospheric  NH3.  Two  different  experiments  were 
performed.  The  first  experiment  involved  the  use  of  the  dual -laser  DIAL  system; 
these  measurements  used  the  less  than  optimal  P(30)  and  P(32)  CO2  laser  lines 
and  LIDAR  returns  from  topographic  targets.  The  second  set  of  experiments  used 
a  single-laser  DIAL  system;  these  measu:  _nts  employed  the  more  optimal  R(30) 
and  R(26)  lines,  but  used  LIDAR  returns  from  a  retroreflector  due  to  the 


McClatchey  et  As  seen  in  Table  I,  significant  absorption  by  NH3 

occurs  in  both  the  10.6  wm  and  9.6  wni  transition  bands  of  the  CO2  laser.  The 
most  suitable  laser  transition  pair  seems  to  be  the  R(30)  transition  at  9,217  pm 
as  the  high-absorption  line  and  the  R(26)  transition  at  9.237  urn  as  the  low- 
absorption  line.  However,  only  one  CO;  laser  could  be  made  to  operate  on 
these  two  lines,  so  that  only  a  single-laser  DIAL  experiment  was  possible  using 
these  "optimal"  lines.  Dual-laser  DIAl  experiments  were  carried  out  using  the 
less  optimal  P(32)  transition  at  10.716  pm  as  the  high-absorption  line  and  the 
P(30)  transition  at  10.693  pm  as  the  low-absorption  line.  These  two  line  pairs 
provide  closely  spaced  frequencies,  which  helps  to  maximize  the  mutual  coher¬ 
ence^^  of  the  two  laser  beams.  It  should  be  noted  that  while  a  dual-laser 
DIAL  measurement  may  be  more  accurate  in  the  determination  of  (Pa/^a*^  ® 

factor  of  as  much  as  2  or  3  over  that  of  a  single-laser  DIAL  measurement, ^ 
this  factor  can  be  more  than  compensated  in  Eq.  (1)  by  the  increased  accuracy 
due  to  the  use  of  an  "optimal"  absorption  line  where  (og*  -  03)  is  large; 
such  is  the  case  for  the  two  line  pairs  for  NH3  detection  presented  in  this 
paper. 

IV.  Ammonia  Absorption  Measurements  and  Calibration 

The  absorption  coefficients  of  NH3  were  measured  using  the  laboratory 
absorption  cell  shown  in  Fig.  1.  The  CO;  lasers  were  tuned  to  the  appropriate 
transitions  and  the  normalized  intensities  of  the  laser  beams  after  passage 
through  the  cell  were  detected  and  averaged  for  1000  pulses.  The  cell  was 
filled  with  a  mixture  of  780  ppm  NH3  in  air.  The  transmittance  values 
obtained  correspond  to  an  absorption  coefficient,  o,  of  60  (atm  cm)“^  and  13 
(atm  cm)“l  for  the  two  high-absorption  R(30)  and  P(32)  transitions, 
respectively.  These  results  are  in  good  agreement  with  those  shown  in  Table  I 
and  the  recent  NH3  absorption  data  of  Hawley,  et  al.^^ 
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th  bP  laser  transitions  available  from  the  CO2  lasers  used.  Secondly,  the 
difference  in  the  absorption  levels  of  the  two  laser  frequencies  used  in  a  DIAL 
measurement  must  be  sufficient  to  yield  a  meaningful  differential  absorption  in 
the  LIOAR  returns  for  accurate  atmospheric  remote-sensing  measurements.  In 
addition,  atmospheric  transmission  and  possible  interference  effects  from  other 
atmospheric  species  must  be  considered. 

The  trade-off  of  these  factors  may  be  seen  in  Eq.  (1),  which  gives  the 
deduced  path-averaged  concentration  of  the  absorbing  molecule,  Ng,  over  the 
range,  R,  as 

Na  =  (Pa/Pa')  +  2  («  -  a')R  /  2(a^‘  -  Og)  R  ,  (1) 

where  Pa  is  the  backscattered  signal  power  normalized  to  the  transmitted 
power,  Og  is  the  molecular  absorption  cross  section,  and  a  is  the  background 
extinction  coefficient  of  the  atmosphere;  the  primed  and  unprimed  parameters 
refer  to  values  at  the  low-  and  high-absorption  frequencies,  respectively. 

As  seen  in  Eq.  (1),  the  accuracy  in  the  measurement  of  Ng  is  increased  when 
(°a'  "  °a)  large,  the  accuracy  in  the  measurement  of  Pa/Pa '  is  high, 
and  (a  -  a')  is  negligible. 

Based  upon  the  above  criteria,  possible  laser  lines  for  the  remote  sensing 
of  NH3  along  with  the  corresponding  absorption  coefficients  of  NH3  and  the 
atmospheric  extinction  are  given  in  Table  I.  The  NH3  absorption  coefficients 
for  these  CO2  laser  lines^^  were  obtained  from  measurements  by  Patty 
et  al.i^^  and  by  Brewer  and  Bruce. The  atmospheric  transmission  and  pos¬ 
sible  interference  from  other  molecules  were  obtained  from  the  AFGL  tapes  of 


[ 
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